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ABSTRACT 

A simplified  Mach  number  scaling  law  is  obtained  for  rotational 
and  broadband  noise  components  of  a model  helicopter  rotor.  The  broadband 
noise  sources  are  further  classified  into  low  frequency  and  high  frequency 
components.  The  scaling  laws  are  based  on  the  geometric  and  performance 
parameters  of  the  rotor  and  characteristics  of  the  flow  field. 

The  existing  theory  of  Lowson  and  Ollerhead  is  used  deriving  the 
conventional  sixth  power  law  for  the  rotational  noise  of  geometrically 
similar  blades  operating  in  similar  flow  environ.  The  knowledge  of  un- 
steady aerodynamics  has  been  exploited  to  yield  analytical  formulation 
for  the  low  frequency  broadband  radiation.  The  ambiguous  state  of  the  art 
regarding  the  origin  and  nature  of  high  frequency  broadband  noise  does  not 
permit  such  a straightforward  scaling  law  for  this  frequency  regime.  The 
vortices  are  assumed  to  be  shed  at  unknown  Strouhal  frequency  and  the 
scaling  law  is  derived  by  simply  integrating  the  blade  sectional  velocity 
over  the  span. 

The  M.I.T.  5X71/2  foot  anechoic  wind  tunnel  facility  was  used 
to  perform  experiments  at  controlled  flow  environ.  The  turbulence  was 
generated  at  the  inlet  of  the  tunnel  with  the  aid  of  biplanar  grids  and 
simultaneous  measurements  of  acoustic  and  turbulence  signals  were  made. 

The  experimentally  obtained  results  are  compared  with  the  computed 
intensities  and  spectra  of  rotational  noise,  low  frequency  broadband 
noise  and  high  frequency  broadband  noise  from  model  rotors. 
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CHAPTER  I 


INTRODUCTION 

Helicopter  rotors  generate  complex  acoustic  signatures.  A variety 
of  mechanisms  are  responsible  for  the  radiated  acoustic  energy  and  only 
a few  of  these  mechanisms  are  understood. 

The  sound  radiated  by  a lifting  rotor  at  low  to  moderate  tip  speeds 
is  essentially  due  to  the  time  varying  pressure  distributions  on  the 
blade.  This,  often  termed  as  "rotational  noise",  was  originally  studied 
by  Gutin  [1],  where  only  steady  azimuthal ly  constant  blade  loads  were 
considered.  These  discrete  frequency  components  occur  at  multiples  of 
blade  passage  frequency  and  the  unsteady  nature  of  the  blade  forces 
prevents  the  harmonics  from  decaying,  as  originally  predicted  by  Gutin. 
Random  noise  formerly  was  called  "vortex  noise",  but  investigators  now 
prefer  "broadband  noise",  since  vortex  shedding  itself  is  not  believed 
to  be  the  principal  mechanism. 

At  higher  tip  speeds  and/or  high  flight  speed,  the  flow  over  the 
blade  section  may  exceed  sonic  velocity  and  result  in  generation  of  a 
local  shock.  These  are  termed  as  "shock  noise"  and  are  characterized  by 
sharp  peaks  in  the  acoustic  pressure  time  history,  resulting  in  highly 
annoying  banging  sounds.  Another  dominant  noise  source  is  due  to  "blade- 
vortex  interaction".  A blade  that  passes  close  to  a wake  vortex  filament 
experiences  [2,3]  a rapid  change  in  the  aerodynamic  loading.  This,  often 
termed  as  "blade  slap",  is  dependent  on  the  geometry  of  the  wake.  Both 
these  noise  sources  are  impulsive  in  nature  and  occur  only  at  certain 
flight  conditions  of  a helicopter.  In  the  absence  of  these  mechanisms, 
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the  discrete  and  broadband  noise  dominate  the  rotor  noise  spectra. 

Whereas  the  "rotational"  or  "discrete"  noise  from  the  rotors  is 
fairly  well  understood,  the  situation  for  broadband  noise  is  most  unsat- 
isfactory. Almost  nothing  can  be  said  unequivocally  about  the  broadband 
noise  and  this  is  simply  because  of  the  difficulty  in  pinpointing  the 
noise  sources  experimentally.  Opinions  regarding  the  origins  and  be- 
havior of  broadband  noise  vary  somewhat  among  investigators,  but  there 
is  general  agreement  that  turbulence  in  the  flow  seen  by  rotor  blades  is 
the  basic  physical  mechanism  responsible  for  broadband  noise. 

Two  similar  bodies  placed  in  an  otherwise  smooth  flow  generate  sound 
most  efficiently  when  one  lies  in  the  turbulent  wake  of  the  other.  The 
first  person  to  record  such  an  effect  appears  to  have  been  G.I.  Taylor 
[4],  who  waved  a toasting  fork  through  the  air  and  observed  a louder 
"singing"  when  the  plane  of  the  prongs  lay  in  the  direction  of  motion 
than  when  it  was  perpendicular  to  the  motion.  A later  experiment  by 
Kramer  [5]  ingeniously  demonstrated  the  Taylor  effect  for  airfoil 
sections  by  rotating  a single  blade  of  high  aspect  ratio  at  zero  angle 
of  attack,  such  that  the  blade  passed  through  its  own  wake.  The  overall 
noise  level  was  about  3dB  higher  with  this  arrangement  than  when  the 
blade  was  given  a certain  angle  of  attack  in  either  direction.  Recent 
experimental  data  for  a full  scale  helicopter  rotor  obtained  by  leverton 
[6]  and  Leverton  and  Pollard  [7]  seems  to  suggest  that  the  broadband 
noise  is  dependent  only  on  the  nature  of  turbulence.  Thus,  the 
turbulence  may  be  due  to  inflow,  the  boundary  layer  of  the  blades,  or 
the  wake  itself  - generating  an  unsteady  lift  on  the  airfoil  and  hence 
the  noise.  Principal  areas  of  uncertainty  include  the  effects  of 


Li 


r 

t 

- 15  - 

velocity  and  the  size  scale  of  turbulence  on  the  intensity  and  spectra 
of  the  broadband  noise.  The  work  of  Leverton  [6]  fails  to  show  any  kind 
of  Strouhal  frequency  scaling  with  velocity  for  full  scale  rotors. 

Much  of  the  recent  impetus  to  study  the  rotational  and  broadbano 
noise  comes  from  improved  data  processing  technique  and  equipment. 
Narrowband  analysis  of  rotor  noise  data  has  shown  discrete  frequency 
components  extending  well  beyond  150  Hz,  which  historically  was  believed 
to  be  the  transition  region  between  rotational  noise  and  broadband 
noise  for  typical  helicopter  rotors.  Work  by  Lowson,  Whatmore,  and 
Whitfield  [10]  and  by  Leverton  [6]  emphasize  that  low  frequency 
broadband  noise  could  have  different  causative  physical  mechanisms 
and  acoustic  behavior  than  high  frequency  broadband  noise.  It  is 
generally  believed  that  the  inflow  turbulence  interaction  with  the 
rotor  blades  is  responsible  for  the  low  frequency  broadband  noise. 

The  high  frequency  component  is  attributed  to  the  vortex  shedding  by 
the  airfoils.  In  Figure  1.1  is  shown  a typical  6%  bandwidth  acoustic 
spectrum  of  a model  rotor.  Identified  in  this  figure  are  the  regions 
of  rotational  noise,  low  frequency  broadband  noise,  and  high  frequency 
broadband  noise. 

To  obtain  the  relationship  between  flow  field  characteristics  and 
the  sound  radiated  lifting  rotors,  it  is  necessary  to  determine  which 
mechanism  is  primarily  responsible  for  the  radiation.  Each  postulated 
mechanism  leads  to  a particular  characterization  of  the  sound  field  and 
of  its  relationship  to  the  flow  field  characteristics.  Empirical 
measurements  and  similarity  arguments  can  then  be  used  to  establish 
a posteriori  verification  of  the  choice  of  the  mechanism.  To  this  end. 
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we  assume  that  (i)  the  steady  and  unsteady  loading  harmonics  are 
responsible  for  radiation  of  rotational  noise,  (ii)  the  intensity  and 
size  scale  of  turbulence  are  the  governing  factors  for  low  frequency 
broadband  noise  radiation,  and  (iii)  the  periodic  vortex  shedding  is 
the  dominant  mechanism  responsible  for  radiation  of  high  frequency 
broadband  noise.  The  major  contribution  of  the  research  effort 
described  in  this  report  is  the  development  of  appropriate  Mach  number 
scaling  laws  for  model  helicopter  rotor  noise  and  experimental 
verification  of  these  laws.  Equally  important  has  been  the  obtainment 
of  a definition  of  the  effects  of  helicopter  performance  parameters 
on  model  helicopter  rotor  noise. 

Here,  an  attempt  has  been  made  to  study  the  influence  of  free- 
stream  disturbances  on  the  acoustic  signatures  generated  by  model 
helicopter  rotors.  Appropriate  theoretical  modelling  for  rotational, 
low  frequency  broadband  and  high  frequency  broadband  noise  are  made 
exploiting  the  existing  knowledge  of  acoustic  radiation  from  a rotating 
dipole  source  and  unsteady  aerodynamics  associated  with  rotating  blades. 
In  Chapter  II  is  presented  a detailed  description  of  the  MIT  Anechoic 
Wind  Tunnel  and  model  rotor  facility,  the  instrumentation  used  for 
acquisition  and  processing  of  various  rotor  noise  components.  Here 
we  have  discussed  the  nature  of  the  artificially  generated  tunnel 
turbulence  and  the  appropriate  theoretical  modelling  of  its  spectrum. 

In  Chapter  III  is  derived  a simple  Mach  number  scaling  procedure  to 
describe  to  rotational  noise  characteristics.  Chapter  IV  contains  a 
formal  theoretical  analysis  of  the  low  frequency  broadband  noise 
radiation.  We  have  compared  the  experimentally  obtained  data  with  the 
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computed  intensity  and  spectra  of  low  frequency  broadband  noise.  In 
Chapter  V the  nature  of  the  high  frequency  broadband  noise  is  discussed 
in  detail.  An  order  of  magnitude  study  has  been  made  to  predict  the 
peak  intensity  of  the  high  frequency  broadband  noise  and  its  location 
in  the  frequency  domain.  The  effects  of  inflow  turbulence  on  the 
high  frequency  broadband  noise  radiation  are  postulated  and  the  results 
are  compared  with  experimentally  obtained  data.  Chapter  VI  provides 
the  scaling  results  appropriate  for  the  radiation  of  low  rotational 
frequency  broadband  noise  and  the  high  frequency  broadband  noise.  We 
have  also  made  an  attempt  to  compare  our  model  results  with  the  full 
scale  results  obtained  by  Leverton  [6].  The  conclusions  of  our 
efforts  are  discussed  in  Chapter  VII. 


CHAPTER  II 


DESCRIPTION  OF  THE  EXPERIMENTS 


2.1  Introduction 

The  M.I.T.  anechoic  wind  tunnel  facility  was  used  to  study  the  ef- 
fects of  controlled  free  stream  turbulence  on  the  low  frequency  and  the 
high  frequency  broadband  noise  radiation  from  a model  helicopter  rotor. 

The  parameters  under  study  were  the  variation  of  helicopter  performance 
characteristics  such  as  number  of  blades,  thrust  and  advance  ratio.  In 
addition,  we  also  studied  the  effects  of  intensity  and  size  scales  of  in- 
jested  turbulence  on  the  intensity  and  spectra  of  broadband  noise.  Tur- 
bulence was  generated  at  the  tunnel  inlet  by  placing  biplanar  grids  of 
varying  size.  In  this  chapter,  the  experimental  apparatus  used  in  ob- 
taining and  analyzing  both  turbulence  and  acoustic  data  are  described. 

The  nature  of  the  injested  turbulence  and  the  appropriate  theoretical 
modelling  of  the  spectrum  for  evaluation  of  the  acoustic  spectra  are 
discussed  at  length. 

2.2  The  M.I.T.  Anechoic  Tunnel  & Rotor  Facility 

The  5x71/2  foot  open  Jet  wind  tunnel  was  used  for  the  study  of 
effects  of  free  stream  turbulence  on  model  rotor  generated  noise.  The 
maximum  speed  of  the  tunnel  is  75  mph  and  the  open  jet  section  of  the 
tunnel  is  enclosed  in  a 12  x 12  x 24  foot  anechoic  chamber.  The  sides  of 
the  chamber  were  covered  with  Cremer  blocks  and  the  floor  of  the  chamber 
was  covered  with  6 inch  thick  polyurethene  foam.  The  anechoic  properties 
of  the  tunnel  were  measured  and  the  acoustic  cut  off  frequency  above  which 
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free  field  conditions  prevail  was  found  to  be  160  HZ.  The  effect  of  the 
shear  layer  of  the  open  jet  on  refraction  and  scattering  of  acoustic  waves 
were  studied  by  using  aeolian  tones  as  sound  source  and  were  found  to  be 
insignificant  under  the  present  test  conditions.  The  details  of  the 
aerodynamic  and  acoustic  calibrations  of  the  wind  tunnel  facility  are 
described  in  References  [11]  and  [12  ]. 

The  model  helicopter  rotor  system  consists  of  a 4.17  foot  diameter 
rotor  and  has  blades  made  of  fiber  glass.  The  rotor  hub  is  designed  to 
take  up  to  8 blades  and  is  connected  to  a thrust  measuring  dynamometer. 

The  dynamometer  consists  of  four  sets  of  BLH  SPB  3-20-35  semi-conductor 
strain  gages  mounted  on  two  sets  of  flexures.  Lebow  slip  ring  was  used 
to  transmit  the  signals  from  the  rotating  dynamometer.  The  details  of 
the  dynamometer  system  and  its  calibration  may  be  obtained  from  Reference 
[12  ].  The  rotor  shaft  is  driven  by  a 5.8  HP  hydraulic  motor  which  has 
a maximum  rpm  of  1200.  The  rotor  shaft  has  a tilt  capability  of  ± 20°. 

The  details  of  the  model  rotor  are  summarized  in  Table  2.1. 

Controlled  turbulence  was  generated  with  the  aid  of  biplanar  grids 
mounted  in  the  upstream  section  of  the  tunnel.  Figure  2.1  shows  a photo- 
graph of  the  model  rotor  system  with  an  upstream  turbulence  generator  in 
the  background.  The  steady  velocity  in  the  test  section  was  set  using 
a pitot  static  probe  and  the  read  out  was  made  on  a vertical  manometer. 

The  pressure  variation  was  minimal  for  the  test  conditions  reported  here. 
The  fluctuating  velocity  signals  were  measured  with  the  aid  of  a hot  wire 
sensor.  In  Figure  2.2  is  shown  a photograph  of  the  set  up.  Identified 
within  the  figure  are  the  location  of  the  probes  for  measuring  steady 
and  fluctuating  velocity  signals. 
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2.3  Instrumentation 

Data  flow  for  all  the  experiments  performed  for  this  study  was  from 
microphones  and  hot  wire  sensors  to  a magnetic  tape  and  later  from  the 
magnetic  tape  to  appropriate  spectrum  analyzers. 

The  acoustic  measurements  were  made  along  a circular  arc  at  a radius 
of  52  inches  from  the  rotor  disc  plane,  perpendicular  to  the  direction 
of  the  wind  and  passing  through  the  axis  of  the  rotor.  Acoustic  signals 
were  measured  using  a 1/2  inch  B & K microphone  type  4133  with  a cathode 
follower  type  2614.  The  signals  were  amplified  by  a B & K microphone 
amplifier  and  were  recorded  on  magnetic  tape  using  a P.I.  4 channel  tape 
recorder  model  6204  in  direct  mode.  The  microphones  were  calibrated 
using  a B & K piston  phone  type  4220.  In  addition  to  the  acoustic  signal, 
a timing  pulse  at  a period  of  1/Bn  was  also  recorded  on  one  of  the  channels 
of  the  tape  recorder.  The  steady  thrust  measurements  were  made  simulta- 
neously with  the  aid  of  a digital  volt  meter.  In  Figure  2.3  is  shown  a 
schematic  of  the  instrumentation  utilized  in  acquisition  of  acoustic  data. 

The  fluctuating  velocity  signals  were  measured  with  an  X-type  hot 
wire  sensor  type  DISA  55P51 . The  probe  was  placed  7.5"  directly  above 
the  rotor  plane  4"  aft  the  axis  of  the  rotor.  Flow  Corporation  constant 
temperature  anemometer  type  900-1  were  used  in  conjunction  with  Flow 
Corporation  linearizers  type  900-4  and  type  4835-C.  An  attenuator  was 
used  with  the  latter  linearizer  in  order  to  adjust  the  sensitivity  of 
both  the  wires  to  be  same.  The  probe  was  calibrated  over  anticipated 
test  velocities  and  the  responses  of  both  wires  were  found  linear.  The 
linearized  signals  were  processed  through  a Flow  Corporation  sum  and 
difference  unit  to  yield  the  longitudinal  and  vertical  components  of 
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velocities.  These  signals  were  further  passed  through  a Flow  Corporation 
suppressor  filter  unit  type  900-3  to  yield  only  the  fluctuating  quantities 
of  the  signal.  All  the  signals  were  constantly  monitored  with  the  aid 
of  a Flow  Corporation  digital  volt  meter  type  900-2D.  The  fluctuating 
signals  were  then  recorded  on  two  F.M.  channels  of  a P.I.  4 channel  tape 
recorder  type  6204.  The  r.m.s.  values  of  the  signals  were  measured  by 
using  HP  3400A  RMS  volt  meters.  In  Figure  2.4  is  shown  a schematic  of 
instrumentation  used  in  acquisition  of  turbulence  data.  The  acoustic 
and  turbulence  measurements  were  made  simultaneously. 

The  measured  acoustic  signals  were  processed  with  the  aid  of  a 
B & K audio  Frequency  analyzer  type  2107  and  a B & K graphic  level 
recorder  type  2305.  Figure  2.5  depicts  a schematic  of  the  instrumenta- 
tion used  for  analyzing  the  acoustic  data.  The  turbulence  data  was  ana- 
lyzed using  a Federal  Scientific  UA-15A  Spectrum  Analyzer  together  with 
a 1015  averager.  In  Figure  2.6  is  shown  a schematic  of  instrumentation 
for  analyzing  the  turbulence  data. 

2.4  Turbulence  Measurements 

The  grids  employed  in  this  study  were  designed  based  on  the  data 

of  Baines  and  Peterson  [13  ].  The  grids  were  biplanar  consisting  of 

bars  of  3/4  inches  with  a mesh  size  of  6 inches  and  bars  of  3 1/2 

inches  with  a mesh  size  of  20  inches.  The  grid  solidity  were  0.23 

and  0.32  respectively.  The  grid  Reynolds  number  based  on  the  lowest 

4 5 

tunnel  velocity  were  9 x 10  and  3 x 10  respectively.  As  noted  in 
Reference  [13  ],  the  turbulence  down  stream  of  a biplanar  grid  becomes 
practically  isotropic  in  an  axial  distance  of  about  5 to  10  mesh  lengths. 
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In  the  present  series  of  experiments,  the  grids  were  located  96  inches 
from  the  axis  of  the  rotor.  This  corresponds  to  ab^'t  16  mesh  lengths 
for  the  smaller  grid  and  4.8  mesh  lengths  for  the  larger  grid.  The 
longitudinal  and  vertical  integral  scales  and  Ag  of  the  grid  generated 
turbulence  were  determined  near  the  axis  of  rotor.  For  convenience,  we 
estimated  from  the  Eularian  Integral  time  scale  t . The  values  of  xg 
were  determined  from  the  extrapolated  zero  intercept  of  the  power  spectra 
of  longitudinal  and  vertical  velocities.  The  length  scales,  then  are 
given  by 
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Figure  2.7  shows  the  typical  variation  of  measured  longitudinal  integral 
scale  in  the  presence  of  as  well  as  in  the  absence  of  turbulence  generating 
grids.  The  longitudinal  scale  is  roughly  independent  of  the  free  stream 
velocity  when  the  grids  are  present.  In  the  absence  of  grids,  the 
measured  Eularian  time  scales  were  quite  large  and  fluctuating.  This 
resulted  in  large  length  scales  that  are  varying  considerably  with  the 
free  stream  velocity.  However,  the  variation  of  the  length  scale  does 
not  seem  to  follow  any  definite  pattern.  For  convenience  in  the 
theoretical  analysis  we  assumed  that  the  scale  of  turbulence  was  in- 
dependent of  the  free  stream  velocity.  To  this  end,  we  obtained  an 
average  value  of  the  scale  which  is  plotted  in  dotted  lines  in  Figure 
2.7.  In  Figure  2.8  is  shown  the  measured  vertical  integral  scale  as 
a function  of  free  stream  velocity.  Once  again,  we  observe  that  the 
fluctuations  in  the  measured  scales  are  large  in  the  absence  of  grids. 

In  Table  2.2  are  given  the  averaged  values  of  longitudinal  and  vertical 
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length  scales  for  the  three  different  cases.  From  the  table,  it  is  seen 
that  the  ratio  of  vertical  scale  to  longitudinal  scale  is  slightly  greater 
than  the  expected  value  of  0.5  for  isotropic,  homogeneous  turbulence. 

In  Figures  2.9  and  2.10  are  the  measured  power  spectral  density  of 
the  longitudinal  and  vertical  components  of  turbulence  as  a function  of 
wave  number  respectively.  In  plotting  these  figures,  Taylor's  hypothesis 
of  frozen  turbulence  has  been  applied  to  convert  the  turbulence  spectrum 
from  an  experimentally  obtained  temporal  frame  to  a spatial  frame.  The 
axial  wave  number  is  related  to  frequency  by  k = 2irf /U  and  this  is 

X o 

nondimensional ized  with  respect  to  k as  3k  /8-rrk  . Here  k is  the  wave 

e x e e 

number  range  of  energy  containing  eddies  and  is  related  to  the  longitudi- 
nal integral  scale  A^  by  kg  = 0.75/A^.  Also  shown  in  the  figures  is 
the  empirical  isotropic  turbulence  spectrum  derived  from  von  Karman's 
interpolation  formula.  The  nondimensional ized  power  spectral  densities 
for  longitudinal  and  vertical  components  of  turbulence  are  given  as 
follows 
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The  measured  data  shows  good  agreement  with  the  von  Karman  spectrum  in- 
cluding a -5/3  power  decrease  at  high  frequencies.  However,  in  the 
absence  of  grids,  the  tunnel  turbulence  seems  to  deviate  slightly  from 


the  predicted  von  Karman  spectrum.  In  Figure  2.11  is  shown  the  variation 
of  fluctuating  axial  and  vertical  components  of  turbulence  as  a function 
of  free  stream  velocity.  From  the  Figure  2.11  we  infer  than  an  average 
turbulence  intensity  of  11%  for  the  larger  grid  and  6.25%  for  the  smaller 
grid  is  very  reasonable.  In  absence  of  grids,  the  tunnel  turbulence 
seems  to  average  to  about  1.7%  of  the  free  stream  velocity. 

2.5  Acoustic  Measurements 

Extensive  measurements  of  the  background  noise  for  the  M.I.T.  anechoic 
wind  tunnel  facility  has  been  made  by  Harris  and  Lee  [ 15  ]. 

However,  it  is  very  likely  that  an  upstream  turbulence  generator  would  alter 
the  background  noise  spectra  of  the  tunnel.  To  ensure  that  the  acoustic 
measurements  were  not  affected  by  such  noise,  a survey  of  the  tunnel  back- 
ground noise  was  made  at  various  tunnel  speeds  with  and  without  the 
turbulence  generating  grids  in  place. 

To  isolate  the  low  frequency  rotational  noise  components  of  the 
rotor  from  the  relatively  high  intensity  tunnel  background  noise,  a periodic 
sampling  technique  was  used.  This  technique  involves  the  use  of  a PAR  wave- 
form eductor.  Model  TDH-9.  The  waveform  eductor  samples  the  repetitive 
input  and  stores  them  in  a 100  channel  capacitance  memory.  After  a suffi- 
cient number  of  sweeps  have  occurred,  the  noise  from  non-repetitive  sources 
(background  noise  in  the  present  case)  will  be  suppressed  since  their  average 
value  will  approach  zero.  The  processed  rotor  data  will  contain  only  the 
rotational  noise.  The  timing  pulse  of  ^ was  used  to  average  the  waveform 
of  the  signal.  The  spectrum  was  obtained  by  using  a Federal  Scientific  UA-15A 
Spectrum  Analyzer  with  a 1015  averager.  In  Figure  2.12,  the  instrumentation 
used  for  processing  the  rotational  noise  data  is  shown.  In  Figure  2.13,  the 
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typical  waveforms  and  the  resulting  spectra  with  and  without  the  use  of 
periodic  sampling  technique  are  shown. 

Since  some  of  our  studies  involve  the  influence  of  free  stream 
disturbances  on  the  low  frequency  and  high  frequency  broadband  noise 
sources,  it  was  necessary  to  ensure  that  (i)  the  rotational  harmonics 
do  not  contribute  significantly  to  the  broadband  noise  generated  by  the 
rotor  and  (ii)  the  background  noise  level  in  the  tunnel  is  low  enough 
to  yield  the  necessary  signal  to  noise  ratio.  The  first  task  was 
accomplished  with  the  use  of  the  periodic  sampling  technique. 

In  Figure  2.14  are  plotted  the  6%  bandwidth  spectra  for  a two 
bladed  rotor  operating  at  672  rpm  and  an  advance  ratio  of  0.3.  The 
3/4"  x 6"  turbulence  generating  grid  was  located  upstream  and  the  Figure 
compares  the  spectra  obtained  with  and  without  waveform  eductor.  The 
rotational  harmonic  obtained  from  the  periodic  sampling  technique  are 
very  distinct  and  they  decay  very  rapidly.  As  can  be  seen  from 
Figure  2.14,  above  200  HZ,  which  we  assume  as  the  transition  from  rota- 
tional noise  to  broadband  noise,  the  levels  of  broadband  noise  are 
significantly  higher  than  those  of  rotational  harmonics. 

The  second  task  was  accomplished  by  measuring  the  acoustic  signals 
at  a given  observer  location  with  and  without  the  rotor  system  operating. 
Figures  2.15  and  2.16  show  the  measured  6%  bandwidth  sound  pressure  level 
for  3/4"  x 6"  grid  and  3 1/2"  x 20"  grid  respectively.  Also  drawn  in 
these  figures  is  the  typical  noise  spectra  with  the  model  rotor  operating 
in  the  tunnel.  From  the  Figures  2.15  and  2.16  it  is  evident  that  the 
background  noise  in  the  wind  tunnel  does  not  affect  the  rotor  generated 
noise. 
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In  the  absence  of  turbulence  generating  grids,  the  background  noise 
of  the  wind  tunnel  seem  comparable  to  the  intensity  of  low  frequency 
broadband  noise.  Hence,  we  have  not  used  this  experimentally  obtained 
low  frequency  broadband  noise  data  in  our  prediction  techniques.  The 
background  noise,  however,  drops  significantly  above  2000  HZ.  Thus 
the  experimentally  obtained  high  frequency  broadband  noise  data  has 
been  used. 


CHAPTER  III 


DISCRETE  NOISE  RADIATION  FROM  MODEL  ROTORS 

3. 1 Introduction 

The  low  frequency  'rotational'  or  'discrete'  noise  radiated  by 
lifting  rotors  is  essentially  due  to  steady  and  unsteady  pressure 
distributions  on  the  blade.  These  periodic  sources  occur  at  integral 
multiples  of  the  blade  passage  frequency  and  the  unsteady  nature  of 
the  forces  prevent  the  higher  harmonics  from  decaying.  The  rotational 
noise  spectra  of  typical  helicopter  rotors  are  strongly  dependent  on 
the  blade  tip  Mach  number,  in  addition  to  other  helicopter  performance 
parameters.  Most  model  rotor  systems  cannot  simulate  the  high  tip  Mach 
numbers  of  real  helicopters.  In  this  chapter,  we  have  attempted  to  use 
the  existing  theories  of  discrete  tone  radiation  from  subsonic  rotors 
[16  , 17,  18]  to  yield  a simplified  Mach  number  scaling  formula  for 
the  intensity  of  rotational  noise  harmonics.  The  contents  of  this 
chapter  are  largely  an  extension  of  recent  work  by  Aravamudan,  Lee, 
and  Harris  QgO-  In  Section  3.2  is  discussed  the  relative  importance 
of  various  source  terms  contributing  to  the  rotational  noise  radiation. 

In  Section  3.3  is  derived  the  theoretical  Mach  number  scaling  formula. 

The  experimental  verification  of  the  scaling  procedure  is  discussed  in 
Section  3.4. 

3.2  Rotation  Noise  Sources 

The  rotational  noise  of  subsonic  rotors  is  mainly  attributed  to 
force  and  thickness  effects.  The  relative  importance  of  force  and  thick- 
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ness  terms  depend  on  the  rotor  blade  Mach  number.  The  source  strength 
of  the  force  term  is  larger  than  that  of  the  thickness  term.  However, 
the  acoustic  efficiency  of  the  thickness  source  radiation  is  far 
greater  than  that  due  to  the  force.  Thus,  which  mechanism  dominates 
the  acoustic  field  is  essentially  a question  of  whether  the  added 
acoustic  efficiency  of  the  thickness  source  can  overcome  its  inherently 
weak  strength.  Since  the  acoustic  efficiency  decays  rapidly  with 
decreasing  Mach  number,  there  is  a rapid  fall-off  in  acoustic 
effectiveness  for  blade  stations  inboard  of  the  tip.  Hence,  in  the 
following  order-of-magni tude  considerations  on  the  force  and  thickness 
effects,  only  tip  regions  need  to  be  considered. 

For  each  source,  its  effective  Mach  number  M = M*  sina/(l  - M ) 

e t or' 

is  maximum  in  a direction  straight  ahead  of  the  rotor.  The  factor 
(1  - M ) is  minimum  for  the  advancing  blade  and  hence  the  sources  in 
the  vicinity  of  the  advancing  blade  contribute  more.  The  front  rotor 
disc  plane  is  chosen  as  the  observer  location  since  the  thickness  noise 
is  most  intense  in  this  location.  We  further  assume  that  the  tip  Mach 
numbers  are  below  the  critical  Mach  numbers  and  hence  compressibility 
corrections  are  omitted.  Hawkings  and  Lowson  [IS]  found  the  acoustic 
source  strength  of  thickness,  direct  force  and  force  derivative  terms 
to  be  of  the  order  of  pg£l2t,  nF^  , and  1/q,  (3Fr/3t).  The  thickness 
source  term  ^ n2t  may  be  rewritten  as  ( P0  Cq  ^ ) (Mt  c/b)(t/c).  Ignoring 
the  directional  properties  of  F , one  can  rewrite  the  second  term  as 
1/2  CfM2(pac9n)  where  M$  = (1  + p)Mt  for  the  advancing  blade.  Based 
on  quasi -steady  approximations,  Hawkings  and  Lowson  [ISO  suggested  an 
order  of  0.1  pQc 0 si  for  force  derivative  term.  For  typical  helicopter 
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rotors,  b/c  ^ 12,  t/c  ^ 0.12,  ^ 0.06,  and  the  order  of  magnitude  of 

these  sources  are  0.01  Mo  r ft,  0.03  M2  p c ft>  and  0-1  P c ft, 

too  S00  00 

respectively.  Based  on  calculations  of  Hawkings  and  Lowson,  the 
approximate  expressions  for  the  acoustic  efficiency  of  each  of  these 
sources  are  given  by  2.6  x 10“  1 ’ 4 / ( 1 - M^3,  2.2  x 102  M®'V(1  - M^2 

and  1.8  x 102  M®‘3/(1  - M^2,  respectively.  Combining  the  source 
strength  and  its  acoustic  efficiency,  the  ratio  of  acoustic  pressure 
due  to  thickness,  force,  and  force  derivative  is  given  as  follows: 

260  M* ’ 1 : 6.6  M2  M2*1  : 18 
e s e 

For  a typical  helicopter  with  an  advance  ratio  of  0.3,  Mg  = 0.6, 

M$  = 0.65,  and  the  second  term  is  of  the  order  of  unity.  The  thickness 
and  force  derivative  terms  are  equally  important  at  a value  of  Mg  = 0.65 
where  both  the  terms  are  of  same  order.  However,  as  the  Mach  number 
is  decreased  to  about  0.60,  the  thickness  term  is  about  4dB  less  than 
the  corresponding  force  derivative  term  and  at  Mg  = 0.5,  it  is  about 
14dB  less.  Hence,  in  our  Mach  number  scaling  formulation,  we  need 
consider  only  force  derivative  terms  and  the  formulation  would  still 
be  valid  for  relatively  high  up  Mach  numbers;  the  tip  Mach  number  is 
limited  by  0.65. 

3.3  Mach  Number  Scaling  of  Rotational  Noise 

Lowson  and  Ollerhead  derived  the  expressions  for  the  far  field 
acoustic  radiation  due  to  a rotating  fluctuating  point  force.  The 
complex  Fourier  coefficient  Cn  of  sound  radiation  from  a source  at  a 
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distance  R from  the  axis  of  rotation  is  given  by: 


Cn=  l 


” n\-(  n-X) 


ill 


xio  4’ 


n!f7  ^xt<j„-x  + <->Xjnn> 


£,ri 


- bXT<J„-X  * l->ijn+X>]  - '4r  [<n->)Jn-X  * 


n+X' J Rr 


-Rff«"-X>Jn-X-<-'>X<"+‘>Vx] 

(3.1) 

In  Equation  (3.1),  the  subscripts  T,  D,  and  C refer  to  Fourier  coeffi- 
cients of  thrust,  drag,  and  radial  forces  corresponding  to  a loading 
harmonic  X and  the  argument  of  Bessel  functions  is  nMy/ri . The 
coordinate  system  is  shown  in  Figure  3.1.  The  effect  of  blade  number 
may  be  included  by  replacing  n by  mB,  where  m is  a positive  integer. 

For  X=0,  we  obtain  the  case  of  steady  blade  loading  and  Equation  (3.1) 
is  reduced  to  the  classical  propeller  noise  study  of  Gutin.  The 
effect  of  forward  velocity  may  be  included  by  replacing  r by  rlO"MQr)* 
It  should  be  noted  that  this  transformation  applied  to  the  retarded 
position  of  the  rotor,  i.e.,  the  position  of  the  rotor  when  the  sound 
was  emitted. 

In  order  to  obtain  the  desired  scaling  formula,  it  is  necessary  to 


nfly 
c r2 

n l 


taxc<j;. 


,bxc(j;. 
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study  the  characteristics  of  the  two  terms  Jn+^  and  Jn_^  appearing 
repeatedly  in  Equation  (3.1).  The  typical  values  of  n and  X are  much 
larger  than  unity  except  for  the  case  of  steady  blade  loading  and  low 
harmonics.  It  has  been  shown  that  the  range  of  loading  harmonic  X 
contributing  significantly  to  the  sound  radiation  lies  between 
n(l  - KM  y/r  ) where  K ^ 1.36.  For  this  range  of  X,  Jn+^  « and 
hence,  the  former  can  be  ignored.  This  result  can  be  explained 
physically  as  follows.  The  Jn+^  term  represents  the  slow  rotational 
mode,  rotating  at  a speed  of  nft/(n+X)  while  the  J . term  represents  a 
mode  rotating  with  a speed  of  nft/(n-X).  Since  the  faster  modes  have 
higher  efficiency,  the  contribution  of  slower  modes  may  be  neglected. 
Further,  the  proportions  of  the  thrust,  drag,  and  radial  forces  of  a 
typical  helicopter  rotor  is  of  the  order  of  10:1:1.  Hence,  over  the 
efficient  region  of  radiation,  the  contribution  of  drag  and  the  radial 
forces  may  be  ignored.  In  fact,  the  numerical  computation  of  the  noise 
radiation  of  each  of  these  sources  given  by  Lowson  and  Ollerhead  [16] 
justify  such  a simplification  except  near  the  plane  of  the  rotor  disc 
where  the  contribution  of  the  lift  component  would  be  zero  and  the 
drag  and  centrifugal  forces  dominate  the  acoustic  field.  Hence,  over 
the  range  of  angles  except  the  plane  of  the  rotor,  with  these 
simplifications,  Equation  (3.1)  reduces  to 

C"  = jo  77  Cia>T  ' > (3-; 

0 1 

For  the  case  of  a helicopter  rotor,  the  observed  sound  is  the  result  of 
the  continuous  distributed  loading  along  the  blade  span.  Therefore, 


I 
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the  power  spectral  density  of  radiated  sound  P(x,n)  at  an  observer 
location  x is  given  by 


P(x,n)  = lim  ^ 
T-*=°  1 


dt 


Cn(x,n)Cn(x,n' )dndn' 


(3.3) 


0 nn' 


where  Cn  is  the  complex  conjugate  of  Cn<  To  this  end,  we  define  the 
following  non-dimensional  parameters.  Let 


5 = n/b 


(3.4) 


Then, 


nMy/r^  = (nf^y/^  )C 


Let, 


^*XT  + ^XT^  = 2 (Vn^2<:  CXT 


Then, 
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AT' 


(3.5) 


where  $ is  a function  to  account  for  the  phase  difference.  Combining 

A 

Equations  (3.3),  (3.4),  and  (3.5),  we  obtain 
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P(x,n) 


r—  mm 

2 11 

nftX 

n 

_ 00  oo  rr 

[y  P b!s?c]V  l 1 

x-o  x -0  g 

c r 

— o » — 

CXTCX'T^';! 


• «,[V*X'  Vx<nV/^>Vx'(nV/r,5,),i5d5 


(3.6) 


i[W] 

The  term  e accounting  for  the  phase  difference  is  independent 

of  the  tip  Mach  number  M^  and  the  observer  location.  The  term  C^.  may 
be  interpreted  as  the  unsteady  sectional  lift  coefficient  and,  in 
general,  is  a function  of  the  blade  geometry.  For  geometrically 
similar  blades,  the  quantities  inside  the  summation  sign  remain 
unchanged  provided  [M  (y/r  )]  = [M.(y/r  )]  . Hence,  the  Mach  number 

t 11  L 12 

scaling  formula  reduces  to 


P(x2,n,Mt  ) = P(x  ,n,Mt  ) 
2 1 
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(3.7) 


where  sina  = x/r  . 

i 

In  the  case  of  forward  flight,  rj  should  be  replaced  by  r:(l  - MQr) 
where  MQr  is  the  component  of  forward  flight  Mach  number  in  the  direction 
of  the  observer.  With  reference  to  Figure  3.1,  we  have: 


Mor  = Mi  cos®  = uMt  cos®  = “ uM^cos^cosa  (3.8) 

Further,  the  factors  x/r*  and  M^y/r^  must  be  replaced  by: 
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With  these,  the  equation  for  power  spectral  density  becomes. 
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(3.10) 


Once  again,  when  the  arguments  of  the  Bessel  function  in  Equation  (3.10) 
are  equal,  the  Mach  number  scaling  formula  becomes: 


P(xa,n,Mt  ) = P(xi,n,M(;  ) 
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A special  case  arises  when  the  observer  location  is  above  the  axis  of 
the  rotor.  In  this  case,  y/r  = 0,  and  hence  the  only  term  in  (3.1) 
contributing  to  the  entire  sound  field  will  be  when  n = X.  Hence,  both 
(3.7)  and  (3.11)  reduce  to: 

P(x2 ,n,Mt  ) = P(xt ,n,Mt  ) 

2 1 

which  is  the  conventional  sixth  power  Mach  number  scaling  law. 

3.4  Experimental  Results  and  Discussion 

To  verify  the  validity  of  the  Mach  number  scaling  formula  derived 
in  the  preceeding  section,  a detailed  experimental  program  was  carried 
out  in  the  MIT  Anechoic  Wind  Tunnel  facility.  The  details  of  the  Wind 
Tunnel  facility  and  the  instrumentation  utilized  in  data  acquisition 
and  data  processing  are  given  in  Chapter  II  of  this  report.  A periodic 
sampling  technique  was  used  in  determining  the  intensity  of  rotational 
harmonics. 

In  Figures  3.2  and  3.3  are  shown  the  results  of  the  Mach  number 
scaling  formula  for  lower  and  higher  rotational  harmonics  of  a two 
bladed  rotor  on  axis.  The  experimental  resuUs  at  a tip  Mach  number 
of  0.15  were  used  to  predict  the  sound  pressure  levels  at  other  Mach 
numbers.  As  can  be  seen  from  Figures  3.2  and  3.3,  the  theoretical 
Mach  number  scaling  formula  given  by  Equation  (3.12)  yields  better 
comparison  with  the  experimentally  obtained  data  for  higher  rotational 
harmonics.  The  discrepancy  of  the  scaling  formula  for  lower  rotational 
harmonics  is  attributed  to  the  fact  that  the  measurements  were  made  in 
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the  acoustic  near  field  because  of  the  physical  limitations  of  the 
wind  tunnel  facility.  In  Figures  3.4  and  3.5  are  shown  the  scaling 
results  for  the  lower  and  higher  harmonics  of  a two  bladed  rotor  off 
axis.  Here  the  experimental  results  of  the  rotor  operating  at  a tip 
Mach  number  of  0.15  and  at  an  angle  of  59°  were  used  to  predict  the 
sound  pressure  levels  at  other  Mach  numbers.  In  Figures  3.6  and  3.7 
are  shown  the  Mach  number  scaling  results  for  a three  bladed  rotor  on 
axis  and  off  axis,  respectively.  Once  again,  we  observe  that  the 
predictions  for  the  lower  harmonics  are  poorly  correlated  as  compared 
to  those  at  higher  rotational  harmonics. 

Figures  3.8  and  3.9  depict  the  directivity  of  lower  and  higher 
harmonics  of  a two  bladed  rotor.  In  Figure  3.10  is  shown  the 
directivity  of  rotational  Karmonics  of  a three  bladed  rotor.  The 
directivity  data  were  obtained  by  positioning  the  microphones  at 
10°  intervals  and  using  the  periodic  sampling  technique  described  in 
Chapter  II.  As  can  be  seen  from  these  figures,  the  directivity  is 
dependent  on  the  harmonic  number  and  becomes  all  the  more  important 
for  the  off  axis  Mach  number  scaling  law  derived  in  Section  3.3. 


CHAPTER  IV 


LOW  FREQUENCY  BROADBAND  NOISE  RADIATION  FROM  MODEL  ROTORS 

4.1  Introduction 

Since  the  broadband  noise  often  appears  to  follow  roughly  a M 6 

0 

dependence,  it  is  natural  to  postulate  that  unsteady,  but  random,  forces 
on  blades  are  responsible  for  its  generation.  Lowson  and  Oiler-head  [16  ], 
Ffowcs  Williams  and  Hawkings  [20  ] and  Morfey  and  Tanna  [21]  have 
theoretically  related  the  radiated  broadband  acoustic  spectrum  from 
a single  blade  to  the  blade  load  spectrum,  without  specifying  how  the 
loads  are  generated.  Their  results  take  no  account  of  blade-to-blade 
correlations  which  might  exist  in  practical  applications.  Lowson  [ 22] 
has  tried  to  account  for  this  by  postulating  a frequency  independent 
spatial  correlation  scale.  Rotational  harmonics  are  generated  by  the 
passage  of  blades  through  steady  asymmetric  velocity  variations  in  the 
plane  of  the  rotor.  The  broadband  noise,  on  the  other  hand,  is  generated 
by  passage  of  blades  through  a turbulent  - and  thus  a randomly  varying 
velocity  field.  The  essential  difference  between  the  harmonic  loading 
and  random  loading  is  the  bandwidth  of  excitation.  The  origin  of  tur- 
bulence, in  most  cases  seems  rather  inconsequential.  It  is  mainly  due 
to  atmospheric  inhomogeneities,  but  upstream  obstructions  such  as  fuse- 
lage and/or  the  wakes  shed  by  the  moving  airfoil  might  very  well  be  the 
culprit.  A lifting  rotor  produces  a mean  downward  velocity  field  which 
draws  the  atmospheric  eddies  through  the  rotor  plane  with  a convecting 
velocity  V . This  random  variation  in  "up  wash"  induces  a random 
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variation  in  the  angle  of  incidence  and  hence  a fluctuating  blade  load. 
As  atmospheric  turbulence  contains  a spectrum  of  wave  number  components, 
the  resulting  loading  spectrum  affects  tha  acoustic  spectrum  over  the 
entire  range  of  frequencies. 

In  this  chapter,  an  attempt  has  been  made  to  study  the  effects  of 
inflow  turbulence  on  the  intensity  and  spectrum  of  low  frequency  broad- 
band noise.  Here,  we  have  exploited  the  existing  theories  of  radiation 
from  subsonic  rotors  and  the  knowledge  of  aerodynamic  transfer  functions 
to  yield  appropriate  formulae  for  the  radiation  of  the  low  frequency 
broadband  noise.  In  Section  4.2  is  described  the  effects  of  blade-to- 
blade  correlation.  In  Section  4.3  is  given  a brief  account  of  the 
formulation.  Section  4.4  is  giver,  an  account  of  blade  loading  due  to 
turbulence.  To  facilitate  comparison  of  experimentally  obtained  results 
with  the  theoretically  predicted  ones,  we  have  chosen  von  Karman's 
spectrum  for  turbulence  upwash.  The  details  of  the  spectrum  and  aero- 
dynamic transfer  functions  are  given  in  Section  4.5  and  the  associated 
numerical  calculations  are  given  in  Section  4.6.  The  numerically  ob- 
tained results  are  compared  with  experimentally  obtained  data  and  a 
detailed  discussion  of  results  are  given  in  Section  4.7. 

4.2  Effects  of  Blade-to-Blade  Correlation 

Many  of  the  possible  mechanisms  for  broadband  noise  generation  by 
a single  stationary  airfoil  in  a moving  stream  have  been  subject  of 
various  investigators.  The  non-uniform  velocity  field  associated  with 
turbulence  leads  to  an  unsteady  blade  load  that  is  dipole  in  nature. 
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Amiet  [23  ] studied  the  resulting  noise  pattern  by  using  the  non-station- 
ary  airfoil  theory  of  Sears  [24  ].  However,  in  the  case  of  rotating 
blades,  the  problem  becomes  more  complex.  The  very  fact  that  the  blades 
are  rotating  can  be  looked  at  in  terms  of  Doppler  shifts  in  frequency 
[16  ].  The  presence  or  absence  of  a blade  at  various  spatial  locations 
may  in  fact  be  interpreted  as  temporal  variations  seen  in  the  observer 
coordinates.  Thus,  a wide  range  of  frequencies  are  generated  at  a 
single  blade  passage  harmonic  mBQ. 

A hovering  or  low  forward  speed  rotor  also  modifies  aqy  ingested 
turbulence  due  to  the  distortion  of  fluid  elements  drawn  into  the  rotor 
plane.  The  eddies,  in  general,  are  elongated  in  the  direction  of  con- 
vection velocity  and  often  are  chopped  several  times  as  they  pass 
through  the  rotor  plane.  For  example,  an  atmospheric  eddy  of  length 
scale  of  about  10  ft.  will  take  the  order  of  one  second  to  be  drawn 
through  a rotor  plane  with  a convection  velocity  of  about  10  ft. /sec. 
Thus,  for  a two  bladed  rotor  at  5 HZ,  about  twenty  blade  passages  will 
encounter  parts  of  the  same  eddy.  Thus  blade-to-blade  correlation  effect 
allows  positive  and  negative  interference  between  acoustic  waves  genera- 
ted by  various  blades.  This  loading  variation  phenomenon  can  also  be 
looked  at  as  a slow  time  modulation  of  the  loading  seen  by  the  blades 
passing  a given  point.  This  low  frequency  amplitude  modulation  will  then 
add  narrow  side  bands  to  each  blade  passing  harmonic.  On  the  other 
hand,  the  smaller  eddies  which  are  associated  with  high  frequency  sound 
spend  less  time  passing  through  the  rotor  plane  because  of  their  smaller 
spatial  extent.  In  the  above  cited  example,  an  eddy  of  length  scale  of 
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a foot  or  less  will  completely  pass  through  the  plane  of  the  rotor  between 
consecutive  blade  passage.  The  resulting  rapid  modulation  of  the  high 
frequency  blade  loading  implies  wider  side  bands  and  hence  the  radiated 
sound  spreads  over  a large  part  of  the  spectrum.  Hanson  [ 25  ] studied 
this  phenomenon  for  the  case  of  a stationary  fan  and  observed  "partially 
coherent'  blade  loading  which  leads  to  partially  coherent  or  narrow  band 
random  noise  that  are  characterized  by  "hay  stacking"  at  higher  harmo- 
nics of  blade  passage.  Homicz  and  George  [26]  derived  expressions  for 
radiated  high  frequency  noise  which  included  blade-to-blade  correlations 
and  inferred  that  the  sums  and  integrals  involved  in  a modal  analysis 
theory  had  to  be  centered  over  a range  of  loading  harmonics,  Doppler  shift 
terms  and  azimuthal  component  terms.  The  large  number  of  terms  required 
a prohibitive  computation  time  for  evaluation  of  the  spectra  beyond 
500  to  1000  HZ  for  typical  Helicopter  rotors.  In  a subsequent  article, 
George  and  Kim  [27]  and  Amiet  [28]  gave  simpler  formulations  to 
specialize  the  radiation  due  to  atmospheric  turbulence.  Our  analysis 
follows  closely  the  one  proposed  by  Homicz  and  George  [26]  and  George 
and  Kim  [ 27  ]. 

4.3  Formulation 

Figure  4.1  depicts  the  rotor  geometry  showing  the  definitions  of 
symbols  and  convection  velocity  of  a turbulent  eddy.  The  time  interval 
At  for  equally  spaced  blades  to  pass  a fixed  point  on  the  rotor  plane 
is  given  by 


i 
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At  = 


2itR„/B 


(4.1) 


On  the  other  hand,  the  time  scale  for  an  eddy  of  wave  length  X^  to  be 
convected  past  a given  point  in  the  rotor  plane  is  given  by 


(4.2) 


Here  the  eddy  wavelength  X ^ is  given  by  2irA/^,  where  £ is  the  turbulent 
wave  number  and  A is  the  turbulence  integral  scale.  Thus  the  ratio  of 
blade  passage  time  to  convection  time  is 


At  _ 1 * vc 
t_  2ir  B Aft 


The  far  field  radiated  sound  for  loading  due  to  a turbulent  upwash  given 


(4.3) 


by  Homicz  and  George  [26  ] is  of  the  form: 


<S(f)> 


GU.n.m)  l l 
l L n m 

Tni  n 


fR  cos 
0 


Jm  M^(Vc/Aft)  - (£  ' nB)  J d? 


(4.4) 


The  two  series  involving  Bessel's  Function  arise  because  of  the  resolu- 
tion of  the  turbulence  in  the  polar  coordinates  and  because  of  the  Doppler 
shifts  in  frequency  due  to  motion  of  the  blades.  The  frequency  f beyond 
which  (4.4)  predicts  a smooth  spectrum  is  given  by 
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f B (1  + Mo/MJ 

o_  ^ c 

S2  2(1  - M^COS  <p  ) (4.5) 

Also,  the  values  of  turbulent  wave  number  which  make  signigicant  con- 
tributions to  the  spectrum  at  a given  frequency  f are  given  by  [26  ]. 

(f/ft)  (1  - M0cos  $)  (f/ft) (1  + M0cos  4) 

< C < 

(Vc/Afi  )(1  + Mo/Mc)  ) (4.6) 

Combining  (4.5)  and  (4.6)  the  lower  limit  for  5 which  makes  significant 
contributions  for  f > f 0 is  obtained  as 

5-  > BAA 

2VC 

Similarly,  for  f » f0,  we  have 

5f  » (4.7) 

*vc 

Using  (4.3),  this  result  may  be  interpreted  as 

fi.  » i f » f 

‘5 

This  result  implies  that  the  small  turbulent  eddies  which  generate 
sound  at  frequencies  f > f 0 will  completely  pass  through  the  rotor  plane 
before  the  next  blade  arrives.  Under  such  circumstances,  we  may  assume 
that  there  is  no  significant  blade-to-blade  correlation  and  hence  the 
total  radiation  is  given  by  simply  adding  the  uncorrelated  sound  power 
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spectral  densities  of  each  of  the  blades.  The  radiation  of  a single 
blade  approximated  by  a rotating  point  dipole  has  been  given  by 
Ffowcs  Williams  and  Hawkings  [20]  as 


<V*,f)>  = 


4P0cor 


l 

n=-oo 


fR.cos  <p 

Dr  (0,  f - nn  )J2  ( °- 


(4.8) 


where  D ( 0 , f ) is  the  power  spectral  density  of  the  dipole  strength  in 
the  radiation  direction. 

For  a dipole  d (0 ) fluctuating  in  amplitude  with  frequency  f and 
rotating  with  a frequency  Q along  a circular  path,  the  force  derivative 
terms  will  certainly  outweigh  the  acceleration  terms  provided  f » J2 
and  the  source  Mach  number  is  subsonic.  Further  discussion  of  this 
point  is  given  by  Morfey  and  Tanna  [21  ].  This  leads  to  2ui f dr (0 ) 
for  a single  frequency  component  force  and  hence  the  intensity  spectrum 
given  by  equation  (4.8) contains  f2.  However,  any  given  component  of 
load  spectrum  radiates  over  a range  of  Doppler  shifted  frequencies. 
Conversely,  the  acoustic  spectrum  at  a frequency  f will  contain  load 
fluctuations  at  frequencies  shifted  from  f by  a range  of  harmonics  Xft. 
The  amplitudes  of  the  Doppler  shifted  contributions  are  seen  to  be 
proportional  to  J2  (fR0cos  <j>/cfl)  which  become  progressively  smaller  for 
n >_  fR0cos  0/c 0 . Thus,  the  primary  Doppler  effects  are  shifts  Af  = nft 
of  less  than  approximately  fRQcos  <t>/ cQ  or  Af/f  < M0cos  $.  Hence  the 
acoustic  spectrum  at  frequency  f is  affected  by  the  dipole  forces  over 
a range  of  frequencies  of  width  2fM0cos<t>  centered  at  f.  A special  case 
of  interest  arrives  when  the  observer  is  directly  above  the  axis  of 
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rotation.  For  this  case,  <t>  = 90°  and  the  radiation  at  any  frequency  is 
affected  by  the  load  fluctuation  only  at  that  frequency  as  in  the  case 
of  a stationary  airfoil.  For  observer  locations  other  than  this,  a 
large  part  of  the  loading  harmonics  contribute  to  the  acoustic  radia- 
ti  on . 

To  this  end,  we  are  left  with  equation  (4.8)  where  the  only  unknown 
is  the  power  spectral  density  of  the  dipole  source  field  D ($,  f-nfi). 

The  dipole  source  field  is  a direct  consequence  of  the  unsteady  lift 
acting  on  the  airfoil  surface.  Hence,  a knowledge  of  the  unsteady 
aerodynamics  is  essential  to  evaluate  the  source  strength  terms  occurring 
in  equation  (4.8).  Further,  we  assume  that  the  rotating  dipole  is  situated 
at  80%  of  the  blade  span,  where  the  maxima  of  the  steady  lift  occurs  [16  ]. 
In  the  following  section,  we  shall  discuss  the  effects  of  the  incident 
turbulence  on  blade  loading. 


4.4  Blade  Loading  Due  to  Incident  Turbulence 

We  next  find  the  unsteady  lift  distribution  on  a blade  moving 
relative  to  a fluid  containing  a homogeneous,  randomly  varying  upwash 
field  w.  We  define  coordinates  fixed  with  respect  to  the  fluid  as 
x1 , y‘,  z1  and  coordinates  with  respect  to  the  airfoil  as  x,  y,  z. 

Figure  4.2  shows  the  chosen  coordinate  system.  Now,  the  four-dimensional 
Fourier  representation  of  the  turbulence  upwash  with  respect  to  the  fluid 
is  given  as  follows 


"(  x1 , t) 


d3k 


F„(k,  f) 


df 


} 

-00 


- 00 


(4.9) 
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where  t is  the  wave  number  vector  of  the  incident  turbulent  field. 
Referring  to  Figure  4.2,  we  see  that  the  blade  moves  with  the  velocity  (J 
with  respect  to  the  fluid  where  (J  = (Ufl,  0,  -Vc).  Hence  the  relation- 
ship between  the  coordinate  systems  x and  x'  is  given  by  x1  = x -Qt. 

At  this  stage,  we  make  Taylor's  hypothesis  that  the  turbulent  upwash 
field  is  frozen  - that  is  it  does  not  vary  during  the  time  of  a blade 
passage.  This  may  be  represented  in  terms  of  the  Fourier  components  of 
the  spectra  as 

Fw(*.f)  - Fw(lc)6(f) 


With  this,  the  turbulent  upwash  field  in  blade  fixed  coordinates  becomes 


w(x,t) 


d3k 


Fw(*)e 


i 2tt[1c  • 


f-Qt] 


The  fluctuation  in  the  upwash  produces  a fluctuation  in  the  lift 
generated  by  the  airfoil.  To  this  end,  we  define  the  elemental  lift 
response  function  of  the  airfoil  as 


AL(t)  = K(kc,k$)  dw(£) 


where  K is  the  aerodynamic  lift  response  function  at  a moving  point  of 
the  rotating  blade  to  a unit  amplitude  of  turbulent  upwash  and  kc  and 
ks  are  the  wave  number  components  in  chordwise  and  spanwise  directions 
respectively.  The  total  lift  generated  by  the  blade  can  then  be  ex- 
pressed as  follows 


L(t) 


d3k  Fw(£)  K(kc,k$)e 


i2ir1c'($t 


Li 


(4.10) 
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Since  the  turbulent  upwash  is  isotropic,  the  coordinate  system  may  be 
rotated  in  such  a way  that  in  the  new  coordinate  system  x is  parallel  to 
Q as  shown  in  Figure  4.2.  This  transforms  equation  (4.12)  to 


(4. 


where  Q = |(J| 


In  order  to  obtain  the  lift  spectral  density,  we  define  a lift  correlation 
function 


RLL(t,t  + t)  = L(t)  L*(t  + t) 


Under  the  assumption  of  a stationary  random  process,  this  becomes  only 
a function  of  the  time  delay  t and  the  power  spectral  density  of  lift 
may  be  expressed  as  follows 


Vf>  4 


! 


RLL(T) 


Substituting  for  L(t)  from  equation  (4.13),  we  obtain 


Vf>  - 
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Noting  the  mathematical  identity  that 


00 

’ 

exp  [i 2tt( kxQ-f )t]dT  = <5(kxQ-f) 

- 00 


we  have  the  one  sided  power  spectral  density  for  f 
00  00 


\L(f)  ■ J 


<lk  dk  W'/Q.  k„,  k ) 


; j 
-00  -00 


» 0 as  follows 
K(f/Q,ky) | 2 


(4.14 


Equation  (4.14)  is  general  in  the  sense  that  any  known  turbulence  upwash 
spectrum  or  an  aerodynamic  transfer  function  may  be  used  to  yield  the 
power  spectral  densities.  To  facilitate  comparison  of  the  theoretical 
calculations  with  experimentally  obtained  results,  we  use  von  Karman's 
turbulence  spectrum.  The  approximate  values  of  the  aerodynamic  transfer 
function  K are  based  on  Osborne's  [ 29  ] extension  of  the  two  dimensional 
Sear's  function. 

The  blade  loading  spectrum  derived  above  is  that  for  a fluctuating 
lift  force  which  is  actually  oriented  at  a slight  angle  to  the  axis 
of  the  rotor.  This  gives  rise  to  a fluctuating  torque  component  in  addi 
tion  to  the  fluctuating  thrust  components.  Since  the  angle  a<-  is  rather 
small  (-6°)  for  typical  helicopter  rotors,  the  contributions  of  the 
torque  components  are  negligible.  Since  the  thrust  is  along  the  axis 
by  definition,  its  contribution  to  force  in  an  observer  direction  is 
simply  Tsin<j>.  The  corresponding  torque  component  is  Tcos^tana^  which 
is  at  least  one  order  less  than  the  lift  component.  Hence  the  dipole 
strength  may  simply  be  chosen  as  that  of  the  lift  component  alone  which 
is  given  by  Dr~<I>LLsin2<j>. 
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The  influence  of  spanwise  correlation  length  of  turbulence  on  radi- 
ated acoustic  intensity  may  be  accounted  for  by  estimating  the  appropriate 
correlation  lengths  for  atmospheric  turbulence  induced  loads.  Variations 
in  the  phase  of  the  spanwise  load  distribution  can  lead  to  significant 
acoustic  interference  effects  for  particular  combination  of  k,  f and  o. 
However,  these  effects  can  be  calculated  only  with  substantially  more 
detailed  knowledge  of  the  spanwise  aerodynamic  loading  than  is  presently 
available.  Further  more,  one  would  expect  that  the  inherent  randomness 
of  the  problem  would  tend  to  mask  such  a behavior. 

4.5  Turbulence  Spectra  & Aerodynamic  Transfer  Function 

Referring  to  the  coordinate  system  described  in  Figure  4.2,  calcu- 
lation of  the  cross  power  spectral  density  of  the  vertical  component  of 
turbulence  w in  the  lateral  direction  is  required.  This  is  a function 
of  the  axial  wave  number  and  spanwise  separation  distance.  The  general 
expression  for  the  energy  spectrum  tensor  yields  for  the  three  wave 
number  spectrum  of  the  vertical  component  as  [30] 


Adopting  the  empirical  von  Karman's  interpolation  formula,  the  three 
dimensional  energy  function  of  the  turbulence,  E(k)  is  given  by 


E(k) 
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where  kg  is  the  wave  number  range  of  energy  containing  eddies.  Integra- 
tion of  (4.15)  over  wave  number  k then  yields  the  two  wave  number  spectrum 
of  the  vertical  component 
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(4.17) 


The  wave  number  range  of  energy  containing  eddies  kg  is  given  by 


k 
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= /r  r(5/6) 
Af  r(i/3) 


0. 7468/Af 


(4.18) 


Substitution  of  (4.18)  into  (4.17)  yields  the  two  dimensional  turbulence 
spectrum  for  the  vertical  component  as 


*ww 


(k  ,k  ) 
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= 0.4548  w^Ai 
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f [1  + 1 .793Af(kJ  + kj)]?/3 


(4.19) 


By  taking  Fourier  transform  of  equation  (4.19)  with  respect  to  k , we 
obtain  a cross  spectrum  at  a given  frequency  which  would  be  measured  by 
two  probes  separated  by  a distance  y in  the  spanwise  direction.  By  in- 
tegrating the  cross  spectrum  over  y,  a spanwise  correlation  length  may 
be  obtained.  Amiet  [23]  has  shown  that  this  correlation  length  for  high 
frequency  limits  may  roughly  be  estimated  as  kx"^.  We  assume  that  each 
such  strip  radiates  independently  of  its  neighbors  with  random  phase 
correlation  between  them.  This  means  that  their  respective  acoustic 
intensities  may  be  added  directly.  The  radiated  intensity  will  then  be 
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proportional  to  the  square  of  the  area,  c2kx~2  and  total  from  bkx  such 
regions  along  the  span  will  scale  as  c2bkx’^.  For  a point  dipole, 
where  kb  «1 , the  appropriate  scaling  factor  is  (cb)2.  Hence  we  see 

X 

that  the  radiation  for  finite  correlation  lengths  is  reduced  by  a factor 
(k  b)’l  With  this,  a spanwise  correction  factor  for  the  estimate  of  the 

A 

acoustic  power  spectral  density,  will  be  (1  + kxb)~^  or  (1  + bf/Q)"^. 

If  we  assume  that  at  every  instant,  each  rotor  blade  behaves  like 
a two  dimensional  airfoil  in  rectilinear  motion  through  a three  dimen- 
sional sinusoidal  upwash  pattern,  the  three  governing  parameters  for 
the  aerodynamic  transfer  function  are  ky,  V and  M0.  Here 


kT  = Ac  + ks  = (k$  + ky^ 

r = tan_1(ky/kx) 


(4.20) 


These  quantities  will  depend  on  the  magnitude  and  direction  of  the  turbu- 
lent wave  vector  £ as  well  as  the  azimuthal  angle  of  the  blade.  During 
each  revolution,  a blade  will  encounter  some  region  in  which  F approaches 
tt/2  or  zero  implying  nearly  steady  but  a three  dimensional  flow  or  nearly 
two-dimensional  but  unsteady  flow  respectively.  In  either  case,  the 
loading  is  reduced  below  that  predicted  by  quasi  steady  two-dimensional 
theory.  Homicz  and  George  [26]  calculated  the  load  response  to  a con- 
vected  sinusoidal  upwash  pattern  for  two  extreme  values  of  F = 0 and 
r=  tt/2  and  inferred  that  the  influence  of  r may  be  ignored  as  compared 
to  the  effect  of  ky  and  MQ.  Osborne  [29]  gave  the  asymtotic  expres- 


sions  for  compressible  two-dimensional  Sears  function  for  r = 0 as 
follows 
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where  S = A - M2  and  |H(k)|  is  the  magnitude  of  the  original  incompres- 
sible Sears  function  [ 24  ] which  to  a very  good  approximation  may  be  given 
as 


1 

I H ( k ) | * (1  + 2Trk)'  7 (4.22) 

Hence  the  aerodynamic  transfer  function  to  be  used  in  conjunction  with 
equation  (4.14)  is  given  by 
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It  was  pointed  by  Mugridge  [31  ],  that  the  net  lift  response  of  a two 
dimensional  thin  airfoil  to  a convected  three  dimensional  upwash  is  less 
than  the  two  dimensional  upwash  response.  If  the  length  scale  of  turbulence 
is  small  compared  to  the  span  of  the  blade,  the  end  effects,  in  general 
are  negligible.* 


* The  authors  wish  to  thank  Professor  Patrick  Leehey  for  pointing  out  the 
inadequacy  of  a two  dimensional  function. 
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For  an  upwash  pattern  given  by 


w exp  [ i ( 2-n-f t - kxx  + k^y)] 


where  k = mb/2  and  kx  = nfc/j,  the  lift  response  function  is  given  by 


clO<x>  ky) 


i k2  + -2  sin  k, 

4bW  >«irn7T7,>  )i] 

x ir  y 7 


(4.24) 


It  may  be  noted  that  as  ky  ->  0,  this  expression  corresponds  to  the  original 
Sears  function.  Order  of  magnitude  estimates  of  the  use  of  Mugridge's 
function  in  place  of  the  two  dimensional  Sears  function  in  equation  (4.23) 
were  made.  This  resulted  in  reduction  of  about  3.5  dB  at  a frequency  of 
200  HZ.  The  difference  between  the  two  responses  above  500  Hz  were  less 
than  a dB. 

4.6  Numerical  Calculations 

Introducing  equations  (4. 23), (4. 19)  in  (4.14)  and  then  in  (4.8), 
we  obtain  the  final  expression  for  the  spectrum  of  low  frequency  broadband 
noise  as  follows 
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f2sin2d> 


2Qc3r2PnO 


+ bf/Q) 


l Dr(f  -nn)J2  ( 


fR.  sin  f 

-A ) 


n =-» 


(4.25) 
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Dr(f) 
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Before  we  proceed  to  the  actual  computation  of  equation  (4.25),  we  need 
to  make  two  approximations.  Firstly,  the  infinite  summation  over  the 
loading  harmonics  has  to  be  reduced  to  a finite  surranation.  George  and 
Kim  [ 27]  have  demonstrated  that  the  resulting  error  would  be  less  than 
one  percent  if  the  series  were  truncated  at 


n . , n 
min  max 


= ± 1-2  J MJcos  4> 


Secondly,  in  the  evaluation  of  the  integral  of  the  unsteady  lift  function 
and  the  upwash  turbulence  spectrum,  the  integration  must  be  truncated  at 
a finite  wave  number  K^.  The  integration  was  performed  numerically  up 
to  a certain  value  Ky  using  Simpson's  rule  and  then  the  upper  limit  of 
integration  was  increased  by  a predetermined  quantity.  The  integration 
scheme  automatically  terminated  the  integration  when  subsequent  increase 
in  the  upper  limit  yielded  less  than  0.1  percent  increase  in  the  value 
of  integral. 


(4.26) 
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4.7  Results  and  Discussion 


Numerical  calculations  were  performed  for  the  measured  values  of 
turbulence  intensities  and  integral  scales  for  various  advance  ratios 
and  tip  Mach  numbers. 

Figures  4.3  and  4.4  compare  the  predicted  and  the  measured  narrow 
band  spectra  for  the  large  and  the  small  turbulence  generating  grids 
respectively.  The  acoustic  spectra  obtained  is  corrected  for  the 
experimentally  obtained  two  bladed  rotor  results  and  with  a 6%  bandwidth 
analysis  by  adding  3 db  to  the  entire  spectrum.  The  experimental  plot 
is  deliberately  terminated  around  3000  Hz,  beyond  which  is  assumed  to 
be  the  region  of  high  frequency  broadband  noise.  As  can  be  seen  from 
Figures  4.3  and  4.4,  the  agreement  between  measured  and  the  predicted 
Spectra  is  better  for  the  larger  grid  with  an  integral  scale  of  about 

5"  than  for  the  smaller  grid  with  an  integral  scale  of  3.3".  The  re- 

sults of  the  rotor  operating  with  no  upstream  turbulence  generator  are 
not  plotted  because  at  these  frequencies  the  background  noise  was  com- 
parable to  the  intensity  of  the  low  frequency  broadband  noise.  In 
Figure  4.5  are  plotted  the  effect  of  advance  ratio  on  the  predicted  peak 
intensity  of  the  low  frequency  broadband  noise.  As  can  be  seen  from  the 
figure,  a doubling  of  advance  ratio  result  in  an  increase  of  peak  sound 

pressure  level  of  about  8 db.  This  result  is  not  surprising  because  of 

the  increased  r.m.s.  value  of  turbulence  in  the  tunnel  with  increasing 
forward  speed.  Since  the  theory  predicts  a 6 db  increase  in  doubling  the 
turbulence  velocity,  only  the  remainder  must  be  thought  of  as  the  true 
influence  of  advance  ratio.  Caution  must  be  exercised  in  comparing  these 
results  with  those  of  full  scale  rotors  since  increase  in  advance  ratio 
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does  not  significantly  change  the  turbulence  encounter. 

In  Figures  4.6,  4.7  and  4.8  are  drawn  the  influence  of  tip  velocity 
on  the  predicted  spectra  of  the  low  frequency  broadband  noise.  The 
corresponding  integral  scales  are  11.8",  5"  and  3.3"  respectively.  An 
examination  of  these  figures  indicates  that  both  the  spectrum  and  in- 
tensity of  low  frequency  broadband  noise  are  strongly  affected  by  the 
scale  of  turbulence  as  well  as  the  rotational  Mach  number. 

In  Figure  4.9,  we  have  plotted  the  peak  intensity  of  the  low 
frequency  broadband  noise  as  a function  of  the  relative  tip  velocity. 

In  plotting  this  figure,  we  removed  the  effect  of  the  increase  in  the 

w2 

r.m.s.  turbulence  levels  by  taking  out  20  log  — from  the  predicted  sound 

W1 

pressure  levels.  The  levels  are  normalized  with  respect  to  the  predicted 
sound  pressure  level  at  a Mach  number  of  0.1.  Also  plotted  in  the  figure 
are  measured  sound  pressure  levels  for  a two  bladed  rotor  operating  with  a 
0.75"  x 6"  grid.  As  can  be  seen  from  the  figure,  a 12  db  increase  per 
doubling  velocity  seems  to  fit  the  predicted  and  the  measured  data  very 
well.  Including  the  effect  of  increase  in  r.m.s.  turbulence  levels  would 
result  in  the  conventional  M®  power  law.  The  location  of  the  peak  in- 
tensity in  the  frequency  domain  increases  with  increase  in  tip  velocities. 

In  Figure  4.10  are  plotted  the  predicted  sound  pressure  levels  as  a function 
of  the  turbulence  integral  scale  for  a fixed  value  of  the  tip  Mach  number 
of  0.25.  Once  again,  we  observe  that  the  location  of  the  peak  intensity 
in  the  frequency  domain  decreases  with  an  increase  in  the  integral  scale 
of  turbulence.  In  Figure  4.11  are  plotted  the  peak  intensity  of  LFBN  as  a 
function  of  the  longitudinal  integral  scale  A^,  of  turbulence.  Here,  the 
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effect  of  the  rotational  velocity  and  the  increase  in  turbulence  intensity 
are  removed.  As  can  be  seen  from  the  figure  the  peak  intensity  is  a very 
weak  function  of  A^.  This  weak  dependence  can  be  readily  explained  by  an 
examination  of  equation  (4.26).  For  relatively  large  values  of  frequency, 
we  may  approximate  the  denominator  of  equation  (4.26)  as 


/(f/Q)2+  k2  ] [1.793  A*  {(f/Q)2+  k2  }] 


Since  the  integral  is  independent  of  A^  we  readily  obtain  a scaling 
factor  based  on  A^  as 


Dr(f)  - Af 


which  would  result  in  a reduction  of  approximately  2 db  per  doubling 
the  integral  scale. 

An  examination  of  Figures  4.3  to  4.11  suggest  that  the  location  of 
the  peak  intensity  of  the  LFBN  in  the  frequency  domain  can  easily  be 
scaled  with  reference  to  a Strouhal  number  defined  as  follows 


S = 

^ \/ 


In  Figure  4.12  is  plotted  the  spectra  of  sound  pressure  levels  for 
various  Mach  numbers  and  size  scale  of  turbulence.  It  is  interesting  to 
see  how  the  entire  data  collapses  to  a single  curve  defining  a normalized 
acoustic  spectra  as  a function  of  Strouhal  number.  The  peak  intensity 
location  seems  to  appear  at  a Strouhal  number  of  about  1.5. 


CHAPTER  V 


HIGH  FREQUENCY  BROADBAND  NOISE  RADIATION  FROM  MODEL  ROTORS 

5.1  Introduction 

In  the  preceeding  chapter  we  discussed  the  possible  mechanisms  of 

low  frequency  broadband  noise  radiation  from  rotating  airfoils  and  arrived 

at  an  acceptable  Strouhal  frequency  scaling  to  predict  intensity  and 

spectrum  of  such  a radiation.  It  was  observed  that  at  high  enough  fre- 
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quencies,  this  radiation  decays  approximately  like  f . Therefore,  in  ab- 
sence of  any  other  mechanism  of  radiation,  the  intensity  of  rotor  generated 
noise  would  be  insignificant  at  frequencies  well  beyond  the  peak  of  low 
frequency  broadband  noise  radiation.  But,  our  experimentally  obtained 
model  data  as  well  as  the  full  scale  rotor  data  of  Leverton  [ 6 ] do  not 
exhibit  such  a decrease  in  the  measured  acoustic  spectra.  In  fact,  the 
evidence  of  existence  of  another  mechanism  which  radiates  efficiently  over 
the  frequency  range  well  beyond  low  frequency  broadband  noise  peak  is 
overwhelming.  In  this  chapter  we  shall  discuss  the  possible  mechanisms 
appropriate  for  such  high  frequency  radiation. 

As  mentioned  before,  the  broadband  part  of  the  rotor  noise  spectrum 
is  attributed  to  the  non-periodic  time  dependent  loads.  These  fluctuations 
in  force  may  arise  as  a result  of  (i)  free  stream  disturbances  generating 
a fluctuating  angle  of  incidence  (ii)  a turbulent  boundary  layer  on  the 
surface  of  the  airfoil  and  (iii)  a spanwise  vorticity  in  the  wake  of  the 
airfoil.  These  mechanisms,  in  general,  interact  through  nonlinear  fluid 
mechanical  processes.  However,  under  certain  conditions,  it  is  possible 
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to  consider  them  as  being  mutually  independent. 

The  interaction  of  the  turbulent  boundary  layer  with  the  trailing 
edge  of  an  airfoil  imparts  momentum  fluctuations  to  the  surrounding  media 
and  the  total  radiated  spectrum  may  be  divided  into  turbulent  boundary 
layer  dominated  components  and  wake  dominated  components.  Chanaud  and 
Hayden  [32  ] have  shown  that  when  the  dimensions  of  the  radiating 
surface  are  small  compared  to  the  acoustic  wave  length,  the  contributions 
of  the  turbulent  boundary  layer  to  the  total  acoustic  intensity  are  usually 
negligible  relative  to  the  contributions  from  the  inflow  turbulence  or 
wake.  In  the  absence  of  signigicant  inflow  turbulence,  the  wake  generated 
noise  would  dominate  the  high  frequency  regime  of  the  acoustic  spectra. 

Several  independent  experiments  by  Hersh  and  Hayden  [33  ],  Paterson, 
Vogt,  Fink  and  Munch  [34  ] and  Clark  [35  ] revealed  that  discrete  tones 
are  emitted  from  isolated  airfoils  under  certain  operating  conditions. 

Over  a wide  range  of  Reynolds  number  based  on  the  chord  of  the  airfoil, 
the  total  radiated  sound  could  be  dominated  by  these  discrete  tones. 
Further,  the  frequencies  of  these  tones  appear  to  be  a well  defined 
function  of  free  stream  velocity  and  form  an  organized  pattern.  There 
is  a considerable  amount  of  difference  of  opinion  regarding  the  Strouhal 
frequency  scaling  of  these  discrete  tones.  Hersh  and  Hayden  [33  ] and 
Paterson,  et  al.  [34  ] have  all  attributed  the  generation  of  these  tones 
to  the  classical  vortex  shedding  along  the  span  of  the  airfoil.  The  shed 
vortices  generate  a periodic  pressure  fluctuation  over  the  airfoil  at  that 
particular  Strouhal  frequency,  resulting  in  a edge  dipole  radiation. 

Tam  [36  ],  on  the  other  hand  believes  that  there  is  very  little  compelling 
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evidence  to  suggest  a direct  relationship  between  the  vorticity  in 
the  wake  and  the  generated  discrete  tones.  He  postulates  that  these 
tones  are  a direct  consequence  of  a self  excited  feedback  loop  of 
aerodynamic  origin  formed  by  the  acoustic  wave,  boundary  layer  and 
the  wake  flow  of  the  airfoil. 

In  our  efforts  to  understand  the  nature  of  High  Frequency  broad- 
band noise,  we  performed  a series  of  experiments  involving  the  influence 
of  performance  parameters  on  the  intensity  and  spectra  of  High  Frequency 
Broadband  Noise  (HFBN).  The  performance  parameters  chosen  were  (i)  the 
steady  thrust  acting  on  the  model  rotor  system,  (ii)  the  advance  ratio 
and  (iii)  the  number  of  blades.  In  addition,  we  also  examined  the 
effect  of  serrations  on  the  nature  of  HFBN.  In  section  5.2  is  described 
the  effects  of  leading  edge  and  pressure  surface  serrations  on  the  in- 
tensity and  spectra  of  the  HFBN.  Section  5.3  describes  the  dependency 
of  HFBN  on  the  chosen  performance  parameters.  In  the  light  of  these 
measurements,  we  discuss  the  Strouhal  frequency  scaling  procedures  adopted 
by  various  investigators  in  section  5.4.  The  appropriate  amplitude  scaling 
formulae  and  the  experimental  results  to  verify  the  proposed  scaling  pro- 
cedure are  discussed  in  section  5.4.  A qualitative  analysis  of  the  influ- 
ence of  the  free  stream  disturbances  on  the  high  frequency  radiation  from 
rotating  airfoil  is  discussed  in  section  5.5.  The  contents  of  section  5.2 
and  5.4  are  largely  due  to  Lee,  Aravamudan,  Bauer  and  Harris  [ 37  ] and 
Aravamudan,  Lee  and  Harris  [ 19  ]. 

5.2  Seration  Studies 


Serrations  made  of  masking  tape  were  used  as  boundary  layer  tripping 
devices.  The  effect  of  serrations  on  the  high  frequency  broadband  noise 


were  studied  by  attaching  serrations  of  various  lengths  and  thicknesses 
along  the  blade  suction  surface,  pressure  surface  and  leading  edge.  In 
all  these  cases  the  rotor  rotational  speed  was  held  constant  at  1 1 OORPM 
and  the  acoustic  measurements  were  made  at  52  inches  (1.32  m)  directly 
above  the  rotor  disc. 

5.2.1  Suction  Surface  Serration 

Single  layer  serrations  of  2 inches  (5cm)  to  15  inches  (38  cm)  in 
length  were  attached  to  the  suction  surface  of  both  the  blades  of  a two 
bladed  rotor.  The  serrations  were  3/8  inch  (9.5  mm)  aft  leading  edge. 
Little  or  no  change  in  the  HFBN  was  observed.  The  thickness  of  the 
tripping  serrations  was  increased  by  increasing  the  number  of  layers. 
Figure  5.1  shows  the  effect  of  increasing  the  thickness  of  serrations 
of  15  inches  in  length  on  the  spectrum  of  HFBN. 

As  can  be  seen  from  the  figure,  little  or  no  change  in  the 
intensity  of  HFBN  occurs  for  serrations  up  to  4 layers  and  significant 
reduction  occurs  for  serrations  of  6 layers.  However,  this  proved  to 
be  of  very  little  practical  interest  since  the  total  thrust  generated 
by  the  rotor  also  decreased.  The  acoustic  spectrum  had  an  even  larger 
change  when  a spoiler  type  serration  of  15  inches  in  length  was  attached 
to  the  blades.  The  HFBN  peak  disappeared  completely  but  a significant 
increase  in  the  low  frequency  broadband  noise  was  observed.  The  thrust 
of  the  blades  also  decreased  sharply  and  the  blades  were  believed  to  be 
fully  stalled.  Paterson  et  al  [34  ] obtained  similar  results  for  the 
case  of  an  isolated  airfoil. 

5.2.2  Leading  Edge  Serration 

The  noise  spectra  of  a two  bladed  rotor  with  serrations  of  various 


lengths  along  the  leading  edge  were  measured.  Figure  5.2  depicts  the 
effect  of  length  of  serration  on  the  HFBN.  The  serrations  were  found 
to  be  effective  in  reducing  the  HFBN,  but  increased  the  low  frequency 
broadband  no'se  considerably.  A 12  inch  (30  cm)  serration  almost 
eliminated  the  HFBN  peak,  but  increased  the  low  frequency  broadband 
noise  by  about  5 db.  This  increase  in  low  frequency  broadband  noise 
may  be  attributed  to  the  increased  level  of  turbulence  and  the  conse- 
quent random  loading  on  the  blades.  A further  increase  in  serrations 
did  not  prove  effective  in  reducing  the  HFBN  peak.  This  seems  to 
suggest  that  the  significant  contributions  to  HFBN  comes  from  the 
outboard  1/2  portion  of  the  blades.  This  finding  is  consistant  with 
that  of  Hersh  and  Hayden  [33  ] for  the  case  of  low  speed  propellers. 

5.2.3  Pressure  Surface  Serration 

Figure  5.3  shows  the  influence  of  serrations  of  various  lengths 
located  3/8"  aft  the  leading  edge  on  the  pressure  surface  of  the  airfoil. 

It  may  be  observed  from  the  figure  that  a 15  inch  serration  reduces  the 
HFBN  by  about  10  db  without  increasing  the  low  frequency  broadband  noise. 
With  a 6 inch  serration  the  peak  level  was  found  to  decrease  by  about 
3 db.  The  surface  pressure  measurements  of  Paterson,  et  al  [ 34]  indicate 
that  the  pressure  surface  boundary  layer  is  laminar  and  the  boundary  layer 
on  the  suction  surface  is  turbulent.  The  serrations  tripped  the  boundary 
layer  on  the  pressure  surface.  Consequently  the  intensity  of  shed  vortices 
were  reduced  and  hence  the  HFBN.  Since  there  was  no  appreciable  increase 
in  the  turbulence  intensity,  the  low  frequency  broadband  noise  remain  un- 
affected. The  boundary  layer  noise  is  increased,  but  since  its  relative 


— ^ 

- 62  - 

magnitude  with  other  sources  is  small,  the  increase  is  not  noticeable. 

Figure  5.4  shows  the  effect  of  a 4 inch  serration  located  at  various 
distances  from  the  tip  of  the  rotor  blade.  As  can  be  seen  from  the  figure, 
the  serrations  are  most  effective  when  they  are  spread  around  75"  of  the 
blade  span.  Location  of  the  serrations  beyond  this  point  seem  to  reduce 
the  efficiency  of  the  process.  Serrations  beyond  12  inches  from  the 
blade  tip  seemed  to  have  no  effect  on  the  HFBN  peak.  This  seems  to  suggest 
that  the  HFBN  peak  comes  from  the  outboard  1/2  blade  span. 

Figure  5.5  shows  the  effects  of  pressure  surface  and  leading  edge 
serrations  on  the  steady  thrust  generated  by  rotor.  As  can  be  seen  from 
Figure  5.5,  the  serrations  had  little  or  no  effect  on  the  total  thrust 
although  it  is  believed  that  the  drag  may  increase  slightly  due  to  the 
existence  of  a turbulent  boundary  layer  on  both  sides  of  the  blades. 

5.3  Parametric  Studies 

< I 

The  effects  of  rotor  thrust,  advance  ratio  and  number  of  blades  on 
the  intensity  and  spectrum  of  HFBN  were  investigated.  The  tip  path  angle 
was  held  constant  at  9°  during  the  entire  set  of  experiments.  Both  narrow- 
band  and  6%  bandwidth  analysis  of  the  data  were  made.  The  directivity  of 
two  and  three  bladed  rotors  were  measured  along  the  rotor  axis  in  a plane 
perpendicular  to  the  wind  direction. 

5.3.1  The  Effects  of  Collective  Pitch  and  Thrust 

Experiments  were  conducted  at  tunnel  speeds  of  44  fps  and  73  fps  with 
a rotor  speed  of  1100  rpm.  The  rotor  collective  pitch  was  set  to  be  0,  2, 

5,  8 and  11.2  degrees  at  the  tip.  Figure  5.6  shows  noise  spectrum  shapes 
for  different  collective  pitches.  A line  representing  the  6th  power  fre- 
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quency  dependence  (f®)  was  drawn  on  each  spectrum.  This  power  law  would 
result  if  each  blade  section  along  the  span  were  radiating  in  a two 
dimensional  manner  according  to  local  speed.  For  low  collective  pitch, 
such  an  approximation  seems  valid.  It  should  be  noted  that  blade-slap 
due  to  blade  vortex  interaction  occurs  at  5 deg.  and  8 deg.  collective 
pitch  settings.  This  implies  that  the  two-dimensional  vortex  shedding 
does  not  prevail  at  these  conditions.  Ham  [38  ] has  shown  that  flow 
separation  is  likely  to  occur  during  blade-vortex  interaction.  The 
separated  flow  may  randomize  the  vortex  shedding  and  consequently  will 
reduce  the  intensity  of  HFBN. 

Figure  5.7  shows  the  6%  bandwidth  analysis  of  HF3N  under  increasing 
thrust  conditions  at  constant  advance  ratio  without  blade  slap.  It  may 
be  noted  from  Figure  5.7  that  neither  the  peak  frequency  location  nor  the 
amplitude  of  HFBN  is  a strong  function  of  steady  state  thrust.  However, 
increase  in  thrust  results  in  a narrowing  of  the  width  of  the  HFBN  hump. 

Also,  a "drag  dipole"  radiation  which  occurs  at  approximately  twice  the 
vortex  shedding  frequency  becomes  more  noticeable  with  increase  of  thrust. 

j 

5.3.2  The  Effects  of  Advance  Ratio 

To  determine  the  effect  of  advance  ratio  on  the  spectrum  and  intensity 
of  HFBN,  the  tunnel  speed  was  varied  from  29.3  ft. /sec.  to  58.7  ft. /sec. 
maintaining  the  rotor  rotational  speed  constant  at  672  rpm.  This  corresponds 
to  a variation  in  advance  ratio  from  0.2  to  0.4  at  an  increment  of  0.05. 

The  mean  thrust  developed  by  the  rotor  was  kept  constant  at  2.84  lbs.  by 
changing  the  collective  pitch  of  the  rotor  blades.  The  tip  path  angle  was 
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maintained  at  9°  by  changing  the  shaft  angle.  Figure  5.8  shows  the  6% 
bandwidth  spectrum  of  the  HFBN  at  various  advance  ratios.  It  may  be  noted 
from  Figure  5.8  that  the  location  of  the  peak  is  not  seriously  affected 
by  increasing  the  advance  ratio.  However,  the  peak  seems  to  broaden  with 
increase  of  advance  ratio  which  would  tend  to  increase  the  spectral 
content  in  the  corresponding  1/3  octave  band.  In  addition,  the  intensity 
of  the  peak  by  itself  is  increasing  with  increasing  advance  ratio. 

In  Figure  5.9  is  shown  variation  of  the  peak  intensity  in  6.3  kHZ 
1/3  octave  band  as  a function  of  increasing  advance  ratio.  A 8 db 
per  doubling  advance  ratio  seems  like  a good  approximation  for  estimating 
the  influence  of  advance  ratio  on  the  intensity  of  HFBN. 

5.3.3  The  Effects  of  Number  of  Blades 

To  determine  the  effect  of  number  of  blades  on  HFBN,  a series  of 
tests  were  conducted  with  2,  3,  4,  5 and  6 blades.  Once  again  the  thrust 
was  maintained  constant  at  3.04  lbs.  by  changing  the  collective  pitch 
of  the  rotor  blades.  The  tip  path  angle  was  maintained  at  9°  by  changing 
the  shaft  angle  and  the  rotational  speed  of  the  rotor  was  held  constant 
at  672  rpm.  Figure  5.10  depicts  the  6%  bandwidth  spectrum  of  HFBN  for 
various  number  of  blades.  Figure  5.10  seems  to  indicate  an  increase  in 
the  intensity  of  HFBN  as  well  as  widening  of  the  spectra  with  increase 
in  number  of  blades.  Since  the  increase  in  number  of  blades  is  likely 
to  cause  increased  vortex  shedding  and  consequently  an  increased  blade- 
shed  vortex  interaction,  these  results  appear  consistent.  Since  the 


increase  in  the  intensity  of  HFBN  is  not  confined  to  a narrow  frequency 
band,  an  octave  band  analysis  was  performed  on  the  same  data.  Figure  5.11 
shows  the  variation  of  the  intensity  of  HFBN  in  6.3  kHZ  octave  band  as  a 
function  of  increasing  number  of  blades.  A 3 db  per  doubling  blade  number 
seems  like  a fair  approximation  to  account  for  the  variation  in  number  of 
blades. 

5.3.4  Directivity  Patterns: 

The  directivity  of  two  and  three  bladed  rotors  were  measured  on  a 
plane  normal  to  the  tunnel  wind  axis  at  10°  intervals.  These  results  are 
plotted  in  Figure  5.12.  Also  shown  in  Figure  5.12  is  the  directivity  of 
a vertical  dipole  source  for  comparison. 

5.4  Scaling  of  High  Frequency  Broadband  Noise 

As  mentioned  ibove,  there  is  considerable  amount  of  difference  of 
opinion  among  various  investigators  regarding  the  Strouhal  frequency 
scaling  of  the  discrete  tones  generated  by  airfoils.  In  this  section, 
we  have  reviewed  some  of  the  existing  scaling  procedures.  In  view  of 
the  experimental  results  obtained  in  Section  5.2  and  Section  5.3,  we 
have  proposed  a simple  scaling  formula  relating  the  frequency  of  discrete 
tones  to  the  geometric  properties  of  the  airfoil. 

The  serration  studies  of  an  isolated  airfoil  by  Hersh  and  Hayden  [ 33 
and  Paterson,  et  al . [ 34]  revealed  that  the  intensity  of  the  discrete 
tones  were  reduced  considerably  by  a boundary  layer  trip  wire  located  at 
about  80"  of  the  chord  on  the  pressure  surface  of  the  airfoil.  From  this, 
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Paterson,  et  al.,  inferred  that  the  pressure  of  a tone  was  associated 
with  an  undisturbed  laminar  boundary  layer  on  the  pressure  side  of  the 
airfoil.  The  frequency  of  shed  vortices  are  then  given  by  [34  ], 
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Using  the  results  of  the  laminar  boundary  layer  for  a flat  plate,  the 
expression  for  the  tone  frequency  was  given  by 
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where  Patterson  et  al  [34  ] chose  the  factor  K as  a parameter  to  fit 
their  experimental  data. 

It  was  noted  by  Tam  [36  ] that  the  observed  growth,  saturation  and 
decay  of  the  intensity  of  the  discrete  tones  was  typical  of  an  acoustically 
excited  feedback  loop  of  aerodynamic  origin.  The  loop  consists  of  an 
unstable  laminar  boundary  layer  on  the  pressure  side  of  the  airfoil  which 
upon  merging  with  the  boundary  layer  in  the  suction  side  forms  the  near 
wake  region  as  shown  in  Figure  5.13.  Forced  by  the  acoustic  waves  that 
impinge  on  the  boundary  layer,  the  disturbances  are  initiated  in  the 
wake  flow  at  the  sharp  trailing  edge  of  the  airfoil  - point  A in  Figure  5.13. 
The  disturbances,  in  the  form  of  large  scale  unstable  waves,  grow  in 
amplitude  as  they  propagate  downstream  along  the  wake.  When  these  dis- 
turbances acquire  sufficiently  large  amplitude,  they  cause  the  wake  to 
vibrate  laterally.  Since  the  acoustic  waves  are  produced  by  the  actual 
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unsteady  lateral  displacement  of  the  wake,  the  noise  source  region  in 
the  present  case  is  highly  localized  - say  point  B in  Figure  5.13.  Thus 
generated  acoustic  waves  propagate  in  all  directions  including  the 
pressure  side  of  the  airfoil  near  the  trailing  edge  where  they  force  the 
boundary  layer  to  oscillate.  In  this  way  the  feedback  loop  is  completed. 

According  to  the  feedback  model  above,  disturbances  are  continuously 
being  initiated  at  the  sharp  trailing  edge  of  the  airfoil.  The  fluid 
dynamic  oscillations  could  be  the  Tollmein-Schl ichting  instability  waves 
within  the  pressure  surface  laminar  boundary  layer.  The  instability 
requirement  of  these  waves  will  now  yield  the  frequencies  of  the  discrete 
tones  generated  by  the  airfoil.  Figure  5.14  depicts  the  stability 
criteria  for  the  Tollmein-Schl ichting  waves  for  a Blasius  profile.  When 
the  boundary  layer  oscillations  are  measured  at  constant  frequency  and 
different  streamwise  locations,  maximum  amplitude  is  found  to  occur  at 
a Reynolds  number  corresponding  to  the  right  branch  of  the  neutral 
stability  curve  shown  in  Figure  5.14.  With  this,  Fink  [39  ] obtained 
the  expression  for  the  frequency  of  tone  radiation  as 


f = (2TT)'1 


(5.3) 


where  6*  is  the  displacement  thickness  and  6<S*/V  is  the  reduced  frequency. 
The  value  of  reduced  frequency  decreases  from  0.17  to  roughly  0.10  as  the 

c £ 

Reynolds  number  based  on  chord  is  increased  from  1 x 10  to  2 x 10  . 

Since  the  Reynolds  number  for  a rotor  blade  varies  along  the  span,  the 
resulting  radiation  is  confined  in  a frequency  band. 
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Hanson  [40  ] made  extensive  measurements  of  the  near  wake  properties 
which  lead  to  the  generation  of  vortex  shedding  sound  from  airfoils.  The 
airfoils  belonged  to  the  NACA  16  series  with  various  trailing  edge  modi- 
fications. He  measured  the  mean  and  fluctuating  base  pressure  and  the 
velocity  profile  in  the  wake  and  related  them  to  the  measured  far  field 
acoustic  spectra.  The  frequency  scaling  of  the  periodic  vortex  noise 
was  given  as  follows 


where 


S*  = 0.0728  (1  - ) 

K ee 


S*  = 1® 
9 us 


R* 

e9 


(5.4) 


Here  SQ*  and  Re^*  are  wake  Strouhal  numbers  and  Reynolds  numbers  based 
on  the  momentum  thickness  9 of  the  near  wake  region.  The  velocity  u$ 
was  given  as 

us  = <5-5> 

where  Cpb  was  the  measured  base  pressure  coefficient. 

5.4.1  Peak  Frequency  Consideration: 


So  far  we  have  considered  various  theories  describing  the  frequency 
scaling  characteristics  of  the  trailing  edge  noise  from  stationary  airfoils. 
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It  should  be  noted  that  no  two  theories  described  herein  predict  the  same 
frequency  and  thus  exists  a controversy  regarding  the  noise  generating 
mechanism.  Several  independent  experimental  investigations  have  confirmed 
the  use  of  a shear  layer  thickness  as  the  appropriate  length  scale  [35,  40, 
45  3 to  predict  the  frequency  of  discrete  tones  from  isolated  airfoils. 

Wake  measurements  for  rotating  airfoils  are  not  available  as  of  now  and 
would  indeed  be  desirable  information.  Lee  [14  ] made  measurements  of 
the  peak  frequency  location  of  the  high  frequency  broadband  no’se  of  the 
present  model  rotor  system,  operating  at  various  speeds.  Based  on  his 
measurements,  he  defined  a Strouhal  number 


where  Vg  was  the  rotational  velocity  measured  at  75%  of  the  blade  span. 
His  experimental  results  lead  to  a Strouhal  number  of  about  1. 

Following  Lee’s  [14]  argument  that  there  is  a need  for  a frequency 
scaling  based  more  on  geometrical  parameters  of  the  system  than  on  the 
fluid  dynamical  characteristics,  we  define  a priori  a frequency  scaling 
law  as 


(5.6) 


The  velocity  V is  chosen  to  be  the  tip  velocity.  Based  on  our  experimental 
results,  the  existence  of  such  a relationship  will  be  demonstrated  later  in 
Section  5.4.3. 
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5.4.2  Pressure  Field  Consideration 

By  developing  the  known  ideas  of  Lighthill  [41  ] concerning  the 
aerodynamic  generation  of  sound  by  unlimited  turbulent  flow,  Curie  [42  ] 
obtained  the  expressions  for  the  acoustic  radiation  from  a surface  S 
immersed  in  turbulent  flow  as  follows: 


P - Po 


— L—  1 — — I TiidT(;, 

4ttCq  Ux^xj  jvolume  r 3xi 


fPj~Pvjvn) 


dxS(y) 


(pv  )SIi 
3t  n r 


(5.7) 


For  an  acoustically  compact  source,  the  far  field  sound  radiatii 
is  given  by  Curies'  extension  of  Lighthill 's  theory. 


X • 

P(x.t)  =-I7c7-7'L  Pi(y»x-x/c0)  ds(y) 
* 0 x2  3t 


(5.8) 


Using  the  coordinate  system  shown  in  Figure  5.15,  we  have  the  instanta- 
neous pressure  due  to  an  element  of  the  airfoil  of  length  AL,  over  which 
the  shed  vortices  are  assumed  to  be  correlated,  given  by 


p ■ - Wc7  C-TT  IF  M (piL> 


The  shed  vortices  cause  a periodic  variation  of  P,.,  the  frequency  being 
given  by  Equation  (5.6).  To  this  end  we  have, 
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Pi  = PqCOS  ( 2irf t ) 


Normalizing  PQ  with  respect  to  the  blade  sectional  velocity  P0v£/2,  we 
have  the  mean  square  pressure  variation  as 


which  may  be  rewritten  in  terms  of  the  advance  ratio  as 


vc  = v [ 1 + y (J)  cos  4>  ] 


(5.10) 


This  has  a maximum  value  at  the  azimuthal  angle  of  $ = 0,  and  substituting 
the  maximum  value  of  (5.10) into  (5.9),  we  obtain 


S2c2p2cos2  9v6  [1+yb/R]4  (AL): 

(AP)2  = 3 

c2  r2  t2 
o 1 


( AP) 2 = S2c2pocos2  9 V2V4(AL)2 


c2  r2  t2 
0 1 


The  blade  sectional  velocity  vc  at  the  chordwise  location  is  given  by 


Vc  = V + UooCOS 


(5.9) 


The  total  pressure  is  therefore  the  sum  of  the  elementary  emitting  surfaces 
or  sources  and  is  given  by 
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SPP^ 


l (Ap)2  = 


sW 

o 

c2t2 


b/AL  v8  (1  + y b/Rk)4(ALk)2 


(5.11) 


Here,  the  subscript  k denotes  the  quantities  corresponding  to  a kth 
element.  Noting  that 

vk  = fiRk  = nkAL 


r 


2 

k 


x2  + R.  2 
k 


cos2ek  = x2/(x2  + Rk2) 


we  can  rewrite  Equation  (5.11)  as  follows: 


S2cV 
SDD(f)  5 2- 

pp  c2t2 

o 


n6(AL)' 


b/4L 

E=1  [1  + ]2 


(5.12) 


Here  we  have  assumed  that  the  correlation  length  AL  is  constant  along 
the  span:  For  the  farfield  solution,  in  general,  Rk/x < 1 ; hence  the 
term  in  the  denominator  of  summation  may  be  approximated  to  unity. 
Further,  since  AL  is  small  compared  to  b,  b/AL  is  a large  number,  say 
N.  Then  the  Equation  (5.12)  becomes 


(f) 


S2c2p2 
o 

c2t2 


ft6  N N N2 

(AL)8  2 k6[l+4p  — + 6y2 — + • • • • ] 

x2  k=l  k k2 


The  summation  is  given  by  standard  formulae  in  [43]  and  for  large  N’s, 
may  be  approximated  as  follows 


r 
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N 


l k6  = N7/7 
k=l 


N 


N7 


N l kS  = -g- 
k=l  0 


I k'  = T 


Hence,  we  obtain  the  total  pressure  field  as 


S (f)  = 
PP 


S2c2p2b  v®AL 
2 t 


C2t2X2 
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r i 4.  1 4 , 42  2 1 
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With  this,  the  Mach  number  scaling  formula  for  the  total  high  frequency 
broadband  noise  generated  by  rotors  in  similar  flow  environ  is  given  by: 


C2V 


VvV  = VvV  * 60  109  ~r^+  20  109  'vt;1 
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and  the  corresponding  peak  frequency  locations  f and  f are  given  by 


1 


(5.13) 


(5.14) 


(5.6).  The  effect  of  blade  number  may  be  included  by  adding  10  log  B /B 

2 1 
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to  Equation  (5.14). 

Mugridge  [44  ] and  Clark,  et  al  [45  ] demonstrated  experimentally 
that  the  correlation  lengths  over  a stationary  airfoil  are  directly  pro- 
portional to  the  displacement  thickness  of  the  turbulent  boundary  layer. 
In  the  case  of  a rotating  airfoil,  the  thickness  of  the  boundary  layer 
is  influenced  by  the  local  Reynolds  number  as  well  as  by  the  spanwise 
flow.  Realistic  estimates  of  correlation  lengths  for  rotating  airfoils 
are  not  presently  available.  However,  since  the  mechanism  is  similar  to 
that  of  a stationary  airfoil,  we  assume  for  scaling  purposes  that  the 
correlation  length  behaves  in  a manner  similar  to  the  boundary  layer 
displacement  thickness  at  a representative  spanwise  location.  Since 
the  displacement  thickness  for  the  turbulent  boundary  layer  varies  like 
c/(Re)0-2.  the  high  frequency  broadband  noise  levels  vary  like  V^5'6. 

5.4.3  Experimental  Results  and  Discussion 

To  verify  the  validity  of  the  frequency  and  amplitude  scaling 
procedures  developed  in  the  previous  sections,  we  performed  a series 
of  experiments  in  the  M.I.T.  anechoic  wind  tunnel  facility.  The  data 
was  analyzed  using  a B & K spectrum  analyzer  with  6%  band  width.  In 
Figure  5.16  are  plotted  typical  HFBN  spectra  for  various  rotational 
speeds  at  constant  advance  ratio  for  a 2 bladed  rotor.  As  noted  above 
in  Section  5.3,  the  advance  ratio  had  no  significant  effect  on  the  peak 
frequency  location,  and  this  justifies  use  of  rotational  tip  velocity 
in  equation  (5.6)  for  Strouhal  scaling  of  shed  vortices. 


- 75  - 


Figure  5.17  shows  a plot  of  peak  frequency  as  a function  of  the  blade 
tip  velocity.  As  can  be  seen  from  Figure  5.17,  the  peak  frequency  is 
directly  proportional  to  the  blade  tip  velocity  with  a proportionality 
constant  of  54.  At  higher  tip  speeds,  the  high  frequency  broadband  noise 
is  spread  over  a narrow  range  of  frequencies  and  hence  there  is  some  scatter 
in  the  measured  data.  Using  maximum  blade  thickness  in  Equation  (5.6),  the 
resulting  Strouhal  frequency  is  1.08.  Although  a wake  characteristic 
length  such  as  displacement  of  momentum  thickness  should  be  the  appropriate 
length  scale  to  be  used,  a simple  relationship  based  on  the  geometrical 
parameter  may  be  very  useful  in  practical  applications. 

In  Figure  5.13,  an  estimate  of  the  peak  sound  pressure  levels  calcu- 
lated using  the  measured  values  of  a two  bladed  rotor  operating  at  550 
rpm  as  a base,  is  plotted.  The  three  bladed  rotor  data  is  normalized 
with  respect  to  the  two  bladed  rotor  data  by  taking  out  lOlogB  /B  from 
measured  sound  pressure  levels.  As  can  be  seen  from  Figure  5.18,  the 
prediction  scheme  presented  in  Section  5.4.2  gives  good  correlations  with 
the  measured  two  blade  and  three  blade  data  except  at  very  low  rotational 
velocities.  It  was  observed  that  at  these  low  velocities,  corresponding 
to  a rotational  speed  of  300  and  400  rpm,  the  measured  high  frequency 
broadband  noise  spectra  showed  intense  peaks.  This  probably  is  due  to 
a fairly  laminar  flow  on  both  sides  of  the  airfoil  resulting  in  more 
intense  shedding  of  vortices.  The  peak  intensity  decreased  slightly  as 
we  increased  the  rotational  velocity  to  about  500  rpm,  and  then  started  in- 
creasing according  to  the  postulated  power  law.  Also  observed  at  these 
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low  velocities  were  multiple  tones  at  twice  and  thrice  the  fundamental 
peak  of  high  frequency  broadband  noise.  These  peaks  were  15  to  20  db 
less  than  the  intensity  of  primary  tone  and  disappeared  as  we  increased 
the  rotational  speed  to  about  500  rpm.  At  present,  we  have  not  been  able 
to  establish  a physical  mechanism  responsible  for  such  secondary  radia- 
tion at  multiples  of  primary  vortex  shedding  frequency. 

When  the  advance  ratio  was  kept  constant,  the  prediction  scheme 
seemed  to  yield  excellent  correlation  with  measured  data.  The  central 
region  in  Figure  5.18  corresponds  to  a variation  in  advance  ratio  which 
again  seems  to  cause  a slight  difference  between  predicted  and  the 
measured  sound  pressure  levels.  Considering  the  number  of  approximations 
made  in  the  scaling  procedure,  the  general  agreement  between  the  theore- 
tical prediction  scheme  and  the  experimental  results  may  be  regarded  as 
highly  favorable. 

5.5  Effects  of  Free  Stream  Turbulence 

It  was  observed  in  Section  5.4  that  the  intensity  of  high  frequency 
broadband  noise  is  dependent  on  the  spanwise  correlation  length  of  shed 
vortices.  It  was  also  assumed  that  these  correlation  lengths  are  propor- 
tional to  the  displacement  thickness  of  the  boundary  layer. 

The  influence  of  free  stream  disturbances  on  the  development  of  a 
turbulent  boundary  layer  on  a flat  plate  has  been  studied  by  McDonald  and 
Kreskovsky  [ 46]  and  Evans  and  Horlock  [ 47].  Using  an  existing  integral 
boundary  layer  calculation  to  predict  the  turbulent  boundary  layers 
developing  in  a turbulent  free  stream,  Evans  and  Horlock  [47  ] concluded 
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that  the  displacement  thickness  decreased  with  increasing  free  stream 
turbulence  intensity.  Extending  their  results  to  the  scaling  procedure 
developed  in  Section  5.4,  we  conclude  that  the  increase  in  the  intensity 
of  free  stream  turbulence  would  in  general  tend  to  decrease  the  intensity 
of  the  high  frequency  broadband  noise.  In  Figure  5.19  is  shown  the 
effect  of  free  stream  turbulence  on  the  6%  bandwidth  high  frequency 
spectrum  of  a two  bladed  rotor  noise.  The  above  prediction  would  result 
in  a decrease  of  5.6  db  for  the  0.75"  x 6"  grid  and  8.1  db  for  the 
3.5"  x 20"  grid  respectively  from  the  measured  intensity  of  high  frequency 
broadband  noise  in  the  absence  of  turbulence  generating  grids.  We  further 
note  that  the  location  of  the  peak  intensity  of  the  high  frequency  broad- 
band noise  remains  unaffected  by  introduction  of  free  stream  disturbances. 
In  Figure  5.20  are  plotted  the  predicted  and  measured  intensities  of  high 
frequency  broadband  noise  as  a function  of  turbulence  intensity.  From 
Figure  5.20  it  is  evident  that  the  influence  of  turbulence  on  the  intensity 
of  HFBN  is  indeed  significant. 

It  should  be  noted  that  the  proposed  scheme  to  include  the  effect 
of  free  stream  turbulence  on  the  intensity  of  HFBN  is  only  qualitative 
in  nature.  The  results  of  Evans  and  Horlock  [ 47]  were  obtained  for  low 
to  moderate  intensities  of  turbulence  (of  the  order  of  4%)  and  the 
validity  of  extending  their  theory  to  the  present  situation  of  high 
intensity  turbulence  is  indeed  argumentative.  But  in  the  absence  of 
accurate  information  regarding  development  of  boundary  layers  on  rotating 
airfoils  in  presence  of  free  stream  turbulence,  only  a crude  approximation 
can  be  made  to  predict  the  nature  of  broadband  noise. 


CHAPTER  VI 


CONCLUSIONS  AND  RECOMMENDATIONS  FOR  FURTHER  RESEARCH 

The  objective  of  this  study  has  been  to  arrive  at  suitable  velocity 
dependent  factors  to  predict  the  various  noise  components  of  a model 
helicopter  rotor.  To  this  end,  we  have  performed  a series  of  experiments 
at  the  M.I.T.  anechoic  wind  tunnel  facility.  The  rotational  and  the 
broadband  noise  components  were  isolated  and  studied  independently.  An 
analysis  based  on  existing  theoretical  modelling  of  unsteady  aerodynamics 
has  been  presented  and  the  computed  results  have  been  compared  with  the 
measured  data.  In  Section  6.1  is  given  the  conclusions  of  our  efforts 
and  Section  6.2  describes  some  of  the  possible  recommendations  for  further 
research  in  this  area. 

6.1  Conclusion 

It  is  postulated  that  the  rotational  noise  from  model  rotors  is  a 
direct  consequence  of  steady  and  unsteady  loading  harmonics  and  the  broad- 
band noise  may  further  be  classified  into  low  frequency  and  high  frequency 
components  because  of  the  inherent  differences  ir.  thier  causative  mechanisms. 
The  low  frequency  broadband  noise  is  attributed  to  the  interaction  of  in- 
gested turbulence  with  the  rotor  blades.  The  high  frequency  broadband  noise 
is  associated  with  the  periodic  vortex  shedding  from  the  rotating  airfoils. 

Analytical  expressions  were  developed  to  determine  the  spectral  content 


of  the  rotational  noise  components  using  existing  theories  of  radiation  from 
rotating  dipole  sources.  An  order  of  magnitude  study  indicated  that  the  con- 
tribution of  force  derivative  terms  outweigh  the  contributions  of  force  and 
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and  thickness  terms  for  tip  Mach  numbers  of  0.6  or  less.  Such  an 
approximation  yielded  simpler  expressions  for  the  complex  Fourier  co- 
efficeints  of  sound  radiation  from  a rotating  point  source.  The  resulting 
Mach  number  scaling  formula  for  rotational  noise  is  as  follows. 
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The  location  in  the  frequency  domain  is  simply  given  by  the  blade  passage 
harmonic,  nBft. 

Theoretical  analysis  of  the  low  frequency  broadband  noise  were  made 
exploiting  the  existing  knowledge  of  unsteady  thin  airfoil  theory  and  the 
characteristics  of  the  inflow  turbulence.  Assumptions  of  isotropy  and 
homogeneity  of  turbulence  were  made  to  make  the  analysis  feasible.  The 
measured  spectrum  of  turbulence  justified  the  validity  of  such  assumptions. 
Further,  it  was  observed  that  the  measured  turbulence  spectra  could  easily 
be  interpreted  in  terms  of  von  Karman's  interpolation  formula.  With  the 
aid  of  this  and  the  two  dimensional  Sears  function,  the  unsteady  lift 
acting  on  the  airfoil  was  computed.  The  resulting  acoustic  intensity  and 
spectra  were  dependent  on  the  flow  and  rotational  Mach  numbers,  intensity 
and  size  scale  of  turbulence.  The  peak  intensity  of  low  frequency  broad- 
band noise  was  found  to  be  varying  with  the  Mach  number  according  to  a 
Mj)  power  law,  and  with  r.m.s  values  of  turbulence  velocity  according  to 
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2 

a w power  law.  The  influence  of  size  scale  of  turbulence  on  the  peak 
intensity  was  found  to  be  insignificant,  approximately  A^"0-33.  However, 
the  location  of  peak  intensity  in  the  frequency  domain  was  found  to  be 
strongly  dependent  on  the  Mach  number  and  the  longitudinal  scale  of 
turbulence.  To  this  end,  we  plotted  the  normalized  acoustic  intensity 
as  a function  of  Strouhal  frequency  fA^/V.  It  was  observed  that  the 
peak  intensity  of  LFBN  occured  at  a Strouhal  number  of  about  1.5.  With 
this,  a scaling  law  can  readily  be  proposed  to  evaluate  the  peak  intensity 
of  low  frequency  broadband  noise.  Under  similar  flow  environ,  the  sound 
pressure  level  of  a rotor  system  2 in  terms  of  a rotor  system  1 will  then 
be  given  as  follows 


SPL 


b c 
2 2 


l 1 


SPL  + 20  log  rr— + 10  log  — 
i b.c. 

i 

[(1  + U2)2  + (Mc  /Mt  )2] 


M, 


+ 40  log  jq-  + 20  log 


Af 

-3.3  log  j-2.  - 20  log 


'0  + Uj)2  + 


2 2 

(Mc  /Mt  H 
1 1 


f 


The  location  of  the  peak  intensity  in  the  frequency  domain  is  given  by 


(6.2) 
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for  respective  rotor  system. 
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To  obtain  a similar  scaling  procedure  for  high  frequency  broadband 
noise  generated  by  a model  rotor,  we  assumed  that  noise  sources  were 
acoustically  compact  and  computed  the  instantaneous  pressure  due  to  an 
element  of  an  airfoil  where  vortices  are  being  shed.  Extending  experi- 
mentally obtained  results  for  the  spanwise  correlation  lengths  for  sta- 
tionary airfoils  to  rotating  airfoils  and  assuming  that  the  correlation 
lengths  vary  like  the  displacement  thickness  of  the  boundary  layer, 
it  was  observed  that  the  peak  intensity  of  high  frequency  broadband  noise 
has  a Vt5'9  factor.  An  expression  which  scales  the  location  of  peak 
intensity  in  the  frequency  domain  was  obtained  based  on  the  rotor  blade 
geometric  parameters.  The  resulting  scaling  laws  for  peak  intensity 
was  found  to  be 
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for  the  respective  rotor  systems.  The  effects  of  intensity  and  size 
scale  of  turbulence  were  less  obvious  in  the  study  of  high  frequency  broad- 
band noise.  One  of  the  effects  of  free  stream  turbulence  is  to  alter  the 
correlation  lengths  of  shed  vortices.  To  this  end,  we  used  an  existing 
integral  boundary  layer  calculation  to  predict  the  turbulent  boundary 
layers  developing  in  a turbulent  free  stream.  The  results  indicated  that 
increase  in  the  intensity  of  free  stream  turbulence  in  general , would 
tend  to  decrease  the  correlation  length,  thus  resulting  in  reducing  the 
intensity  of  high  frequency  broadband  noise. 

The  scaling  laws  obtained  for  the  intensities  of  rotational  noise 
components  were  experimentally  verified  over  a range  of  operating  conditions. 
The  measured  data  showed  general  agreement  with  the  scaling  procedure  except 
at  low  rotational  harmonies.  It  is  believed  that  the  geometric  and  acoustic 
near  field  effects  might  have  caused  some  ambiguity  in  the  measured  rotor 
noise  spectra  at  very  low  frequencies.  The  low  frequency  broadband  noise 
prediction  procedures  were  experimentally  verified  under  various  operating 
conditions  of  the  rotor  and  varying  intensity  and  size  scale  of  turbulence. 

A comparison  of  predicted  peak  frequency  and  sound  pressure  levels  with 
experimental  data  showed  good  agreement  except  for  the  case  of  effect  of 
turbulence  on  the  intensity  of  high  frequency  broadband  noise.  Ir.  this 
case,  the  data  showed  general  agreement  with  the  prediction  technique,  but 
the  measured  intensity  levels  were  higher  than  the  predicted  levels.  We 
further  add  that  improvements  in  modelling  more  of  the  fluid  mechanical 
aspects  of  the  flow  field  are  desired  to  place  this  scaling  procedure  on 
a firmer  basis. 
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With  the  knowledge  of  the  variation  of  the  flow  parameters  on  the 
intensity  and  spectrum  of  low  frequency  and  high  frequency  broadband 
noise,  it  is  indeed  tempting  to  compare  the  scaling  procedure  with  ex- 
perimentally obtained  full  scale  rotor  results  [ 6 ].  Since  the  r.m.s. 
turbulence  velocity  and  size  scale  are  unspecified,  we  estimated  their 
values  based  on  Leverton's  [ 6 ] data  for  a full  scale  rotor  operating 
at  160  rpm  and  generating  a thrust  of  3050  lbs.  The  peak  intensity  was 
measured  to  be  72  db  and  the  frequency  was  250  HZ.  With  this,  we  esti- 
mated the  integral  scale  of  turbulence  as  2.8  ft.  and  the  r.m.s.  turbu- 
lence velocity  as  1.52  ft. /sec.  These  values  were  used  to  predict  the 
intensity  of  low  frequency  broadband  noise  for  the  same  rotor  operating 
at  260  rpm  generating  a thrust  of  700  lbs.  The  model  rotor  used  in  both 
scaling  procedures  was  operating  at  750  rpm,  generating  a thrust  of  3.8 
lbs.  and  the  r.m.s.  velocity  and  size  scale  of  turbulence  were  4.1  ft. /sec. 
and  0.42  ft.  respectively.  This  resulted  in  a prediction  of  80.4  db  at 
a frequency  of  400  HZ.  The  measured  values  were  found  to  be  77  db  at  a 
frequency  of  350  HZ. 

Caution  must  be  exercised  in  comparing  the  high  frequency  broadband 
noise  results  of  the  model  rotor  system  with  those  of  full  scale  tests 
because  of  vast  differences  in  Reynolds  number.  At  high  enough  Reynolds 
number,  there  is  very  little  evidence  that  full  rotors  sustain  a laminar 
boundary  layer.  To  determine  the  relevance  of  the  present  study  to  actual 
rotor  noise,  we  computed  the  intensity  and  spectrum  of  the  full  scale 
rotor  of  Leverton  using  our  model  rotor  running  at  1000  rpm.  Figure  6.1 
shows  the  predicted  and  measured  frequencies  for  the  high  frequency 
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broadband  noise  peak.  In  Figure  6.2,  the  estimated  sound  pressure  levels 
are  compared  with  the  measured  Leverton's  data  after  accounting  for  atmos- 
pheric absorption.  As  can  be  seen  from  Figures  6.1  and  6.2,  the  high 
frequency  broadband  noise  scaling  procedure  predicts  the  location  of  the 
peak  in  the  frequency  domain  accurately,  but  fails  to  predict  the  proper 
sound  pressure  levels.  Since  the  estimated  turbulence  intensity  is  very 
low,  we  attribute  this  discrepancy  to  the  vast  difference  in  the  Reynolds 
number.  The  Reynolds  number  of  the  full  scale  rotor  was  6.62  x 10®  as 

c 

opposed  to  the  Reynolds  number  of  our  model  rotor  of  2.3  X 10  . The  tail 
rotor  data  of  Paterson  et  al  [34]  corresponds  to  a Reynolds  number  of 
2.5  X 10®.  A Sikorsky  5 bladed  rotor  operating  at  a design  rpm  of  1243 
on  a rotor  whirl  test  stand  radiated  90  db  at  10  kHZ.  Using  the  proposed 
prediction  scheme  and  using  the  data  of  our  model  rotor  operating  at 
550  rpm,  we  were  able  to  obtain  a sound  pressure  level  of  90  db  at 
9800  HZ. 

The  analysis  of  Homicz  and  George  [26]  has  placed  the  blade-turbulence 
interaction  noise  modelling  on  a firmer  theoretical  basis  than  existed 
previously  for  the  generation  of  broadband  noise.  But,  their  claim  that 
their  theory  can  be  used  to  predict  the  high  frequency  broadband  noise  is 
unrealistic.  Our  analysis,  being  a direct  consequence  of  the  method  of 
Homicz  and  George  yields  good  comparison  with  the  measured  model  rotor 
data  for  low  frequency  broadband  noise.  But,  their  analysis  predicts  a 
constant  decay  of  the  acoustic  field  with  increasing  frequency,  which  does 
not  happen  in  either  model  or  full  scale  rotor  systems.  We  have,  therefore, 
used  the  classical  vortex  shedding  theory  in  conjunction  with  the  blade- 
turbulence  interaction  theory  to  predict  the  total  acoustic  spectrum  of 
rotor  noise.  The  principal  area  of  uncertainty  lies  in  the  inadequacy 
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of  classical  vortex  shedding  theory  to  predict  the  noise  generated  at 
high  Reynolds  number.  A better  understanding  of  the  fluid  dynamics  of 
the  rotating  blades  are  necessary  to  yield  more  realistic  prediction 
tehcniques  [19]. 

Our  efforts  in  modelling  the  rotational  noise  components  of  the 
model  rotor  and  understanding  the  effects  of  free  stream  turbulence  on 
the  nature  of  high  frequency  and  low  frequency  broadband  noise  components 
have  been  greatly  benefitted  from  the  detailed  experimental  task  performed 
by  Lee  [14].  The  utility  of  our  results  for  scaling  of  low  frequency  and 
high  frequency  broadband  noise  remains  as  demonstrated  in  the  above  com- 
parison of  predictions  and  data.  The  scaling  procedures  developed  in  this 
study  may  be  considered  as  a first  order  attempt  in  responding  to  the 
challenge  of  Mach  number  and  Reynolds  number  scaling  of  full  scale  rotors. 

6.2  Recommendations  for  Further  Research: 

The  objective  of  every  experimental  researcher  is  to  compile  data  at 
various  laboratory  operating  modes  and  to  postulate  prediction  schemes 
applicable  to  full  scale  operating  conditions.  In  this  report,  we  have 
presented  the  results  of  the  detailed  study  of  the  acoustic  characteristics 
of  the  M.I.T.  4.17  foot  diameter  model  rotor  system.  The  tip  Mach  number 
and  Reynolds  number  (based  on  chord)  were  limited  to  0.25  and  2.5  X 105 
respectively  because  of  the  limitations  of  the  rotor  operating  system. 

The  next  logical  step  will  be  to  explore  the  validity  of  the  present  pre- 
diction procedures  at  higher  tip  Mach  numbers  and  Reynolds  numbers.  The 
principal  areas  of  uncertainties  in  extending  the  present  analysis  are  as 
follows : 
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(i)  The  relative  importance  of  force  and  thickness  terms  in  the 
radiation  of  rotational  noise  components  must  be  reconsidered.  The 
order  of  magnitude  analysis  presented  in  Chapter  III  of  this  report 
indicated  that  the  thickness  and  the  force  terms  are  comparable  at 
tip  Mach  numbers  of  about  0.65.  Recently  Farassat  and  Brown  [48], 
using  a time  domain  approach  showed  that  the  thickness  noise  dominated 
the  acoustic  time  history  of  a rotor  operating  at  an  advancing  blade 
Mach  number  of  0.92. 

(ii)  The  effects  of  free  stream  turbulence  on  the  intensity  and 
bandwidth  of  rotational  harmonics  need  to  be  studied  on  a formal  basis. 
It  is  generally  accepted  that  the  turbulence  increases  the  blade  to 
blade  correlation  [26]  and  add  side  bands  at  blade  passage  tones. 

This  effect  was  more  pronounced  for  the  case  of  a stationary  fan 
studied  by  Hanson  [25]. 

(iii)  The  ambiguous  state  of  the  art  regarding  the  frequency  and 
amplitude  scaling  of  high  frequency  broadband  noise  needs  to  be 
clarified.  The  aerodynamic  sources  of  high  frequency  broadband  noise 
are  strongly  dependent  on  the  operating  Reynolds  number  and  has  been 
subject  of  various  investigations  [9,  33,  34,  35,  36].  More  experi- 
mental data  over  a wide  range  of  Reynolds  numbers  is  needed  to  place 
the  theory  of  high  frequency  broadband  noise  on  a firmer  basis. 

It  is  believed  that  blade  surface  pressure  measurements  would  greatly 
enhance  our  understanding  of  the  aerodynamic  noise  generating  mechanisms. 

The  simultaneous  measurements  of  aerodynamic  and  acoustic  signals  will 
facilitate  proper  cross  correlation  analysis. 


r 
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There  is  very  limited  experimental  data  available  for  determining 
the  impulsive  noise  characteristics  of  a helicopter  rotor.  The  blade- 
vortex  interaction  noise,  often  termed  as  "blade  slap"  and  the  shock- 
interaction  noise  often  termed  as  "blade  pop"  are  two  major  sources  of 
impulsive  noise.  The  blade  slap  occur  only  at  certain  operating  conditions 
such  as  a decent,  or  a low  advance  ratio  flight.  Experiments  are  being 
conducted  at  the  M.I.T.  anechoic  tunnel  facility  and  by  Schmitz  at  the 
NASA  Ames  Directorate  [49],  to  gain  more  insight  to  the  characteristics 
of  glade  vortex  interaction  noise.  The  shock  noise  or  "blade  pop"  is 
often  due  to  presence  of  unsteady  shock  waves  on  the  blade  surface.  They 
are  highly  directional,  loud  and  high  pitched  banging  sound  occuring  at 
blade  passing  frequency.  Recent  theoretical  study  by  Kitaplioglu  and 
George  [50]  indicate  that  the  shocks  under  proper  circumstances,  are 
likely  to  make  an  important  contribution  to  helicopter  noise.  We  believe 
that  acquiring  experimental  data  on  shock  noise  would  be  necessary  to 
suppliment  the  theoretical  study  already  in  progress. 
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Table  2.1 

MODEL  ROTOR  CHARACTERISTICS 


Radius  (b) 

25  inches  (64.1  cm) 

Chord  (c) 

2 inches  (5.13  cm) 

Number  of  blades  (B) 

1 to  6 

Section 

NACA  0012 

Twist 

- 8° 

Shaft  tilt  Capability 

±20° 

Maximum  RPM 

1200 

Testing  RPM 

Variable,  from  400  RPM  to  1100  RPM 

Lead-lag 

None 

Cyclic  Pitch 

None 

Collective  Pitch 

By  adjusting  pitch  of  individual  blade 
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Table  2.2 


CHARACTERISTICS  OF  TURBULENCE 


No  grid  3.5"  x 20  " grid  0.75" 
(8.97  cm  x 51.28cm)  (1.92cm 


Longitudinal  Scale  - inches  (cm) 
Vertical  Scale  - inches  (cm) 


Ratio  of  A 


g/Af 


11.8(30.26) 
7 (17.95) 
0.59 

1.7% 


5 (12.82) 
3.2  (8.21) 
0.64 

11% 


3.3 

2.0 

0.61 


x 6"  grid 
x 15.38  cm) 

(8.46) 

(5.13) 


Intensity 


6.25% 
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6%  BANDWIDTH  SPECTRUM  OF  MODEL  ROTOR  NOISE 
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PHOTOGRAPH  OF  THE  ROTOR  ASSEMBLY  WITH  UPSTREAM  TURBULENCE  GENERATOR 


PHOTOGRAPH  OF  THE  SETUP  FOR  MEASURING  ACOUSTIC  AND  TURBULENCE  DATA 


Flow  Corporation 
Linearizers 


SCHEMATIC  OF  INSTRUMENTATION  FOR  ACQUIRING  TURBULENCE  DATA 
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FIGURE  2.5 

SCHEMATIC  OF  INSTRUMENTATION  FOR  PROCESSING 
ACOUSTIC  DATA. 
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X-Y  PLOTTER 


FIGURE  2.6 

SCHEMATIC  OF  INSTRUMENTATION  FOR  PROCESSING 
TURBULENCE  DATA 


EFFECT  OF  FREE  STREAM  VELOCITY  ON  THE  INTEGRAL  SCALE  OF  THE  AXIAL  COMPONENT  OF  TURBULENCE, 
▼ NO  GRID,  # 3.5"  X 20"  GRID,  A 0.75"  X 6"  GRID. 


EFFECT  OF  FREE  STREAM  VELOCITY  ON  THE  INTEGRAL  SCALE  OF  THE  VERTICAL  COMPONENT  OF 
TURBULENCE,  ¥ NO  GRID,  • 3.5"  X 20"  GRID,  ■ 0.75"  X 6"  GRID. 


NORMALIZED  POWER  SPECTRUM,  4 <£uu  (kx  )/<£uu(0) 
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FIGURE  2.9 


ONE  DIMENSIONAL  SPECTRUM  OF  THE  AXIAL  COMPONENT  OF  P°ID  TURBULENCE 


(kx  ) /<f}^yf.O) 


ONE  DIMENSIONAL  SPECTRUM  OF  THE  VERTICAL  COMPONENT  OF  GRID  TURBULENCE 


CITY,  FT/SEC 
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X-Y  PLOTTER, 


FIGURE  2.12 


SCHEMATIC  OF  INSTRUMENTATION  FOR  PROCESSING  ROTATIONAL  NOISE. 


FIGURE  2.13  - ILLUSTRATION  OF  EFFECT  OF  PERIODIC  SAMPLING  TECHNIQUE  - (a)  WAVEFORM  OF  RAW  DATA 
(b)  CORRESPONDING  NARROWBAND  SPECTRA  FROM  0 - 200  Hz  (c)  PROCESSED  WAVEFORM  and 
(d)  SUBESQUENT  NARROWBAND  SPECTRA 


COMPARISON  OF  ROTATIONAL  NOISE  AND  LOW  FREQUENCY  BROADBAND  NOISE  AT  HIGH 
FREQUENCIES  FOR  .75"  X 6"  TURBULENCE  GRID,  B = 2,  0 = 672  RPM,  U = 44  FPS. 


6%  BANDWIDTH  BACKGROUND  SPECTRA  AT  VARIOUS  TUNNEL  SPEEDS  FOR  0.75"  X 6" 
TURBULENCE  GRID,  TOP  CURVE;  ROTOR  GENERATED  NOISE  FOR  B = 3,  fi  = 1000  RPM 
and  U = 51  FPS. 
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6 % BANDWIDTH  BACKGROUND  SPECTRA  AT  VARIOUS  TUNNEL  SPEEDS  FOR  3.5"  X 20" 
TURBULENCE  GRID,  TOP  CURVE;  ROTOR  GENERATED  NOISE  FOR  B = 3,  R = 1000  RPM 
and  u = 51  FPS. 
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SPL  db  (ref  0.0002  n bars ) 
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FIGURE  3.4 


MACH  NUMBER  SCALING  FOR  LOWER  HARMONICS 
OF  A TWO-BLADED  ROTOR  OFF  AXIS 


SPL  db 
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HARMONIC  NUMBER 

FIGURE  3.5 

MACH  NUMBER  SCALING  FOR  HIGHER  HARMONICS 
OF  A TWO-BLADED  ROTOR  OFF  AXIS 
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SOUND  PRESSURE  LEVEL  db  (ref  0.0002  M bars) 
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FIGURE  3.6 

MACH  NUMBER  SCALING  FOR  ROTATIONAL  HARMONICS 
OF  A THREE-BLADED  ROTOR  ON  AXIS 


SOUND  PRESSURE  LEVEL  db  (0.0002  h bars) 


FIGURE  .3.2 

MACH  NUMBER  SCALING  FOR  LOWER  HARMONICS 
OF  A THREE  BLADED  ROTOR  OFF  AXIS 


FIRST  HARMONIC  U =64.5  Ft/sec 

■ SECOND  HARMONIC  £1  = 672  RPM 

■ THIRD  HARMONIC 

• FOURTH  HARMONIC 
. FIFTH  HARMONIC 
- SIXTH  HARMONIC 


FIGURE  3.8 

rRKTW ITTUTTOWlOARMON I CS  “OF  A'TWO-  BLAO ED- ROTOR- 


DIRECTIVITY  OF  HIGHER  HARMONICS  OF  A TWO-BLADED  ROTOR 


EFFECT  OF  MACH  NUMBER  ON  PREDICTED  AND  MEASURED  LOW  FREQUENCY  BROADBAND 
NOISE  FOR  THE  LARGER  GRID  WITH  A A-  = 5.0  INCHES 


EFFECT  OF  MACH  NUMBER  ON  PREDICTED  AND  MEASURED  LOW  FREQUENCY  BROADBAND 
NOISE  FOR  THE  SMALLER  GRID  WITH  A Af  = 3.3  INCHES 


PEAK  SPL  IN  dB,  ref.  0.0002  Mbars 


| 
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FIGURE  4.5 

EFFECT  OF  ADVANCE  RATIO  ON  THE  PEAK  INTENSITY 
OF  LOW  FREQUENCY  BROADBAND  NOISE 

_ _ _ _ J 


FIGURE  4.6 

EFFECT  OF  MACH  NUMBER  ON  PREDICTED  LOW  FREQUENCY  BROADBAND  NOISE  SPECTRA  FOR 
A.  = 11.8  INCHES 


EFFECT  OF  MACH  NUMBER  ON  PREDICTED  LOW  FREQUENCY  BROADBAND  NOISE  SPECTRA 
FOR  A*  = 3.3  INCHES 


SOUND  PRESSURE  LEVE 


Mt  LINE 


PREDICTION: 

■ Af  = 3.3" 

• Af  = 5.0" 

0 Af  =11.8" 

1 EXPERIMENTS 


EFFECT  OF  TURBULENCE  SCALE  ON  THE  PEAK  INTENSIGY  OF  LOW  FREQUENCY  BROADBAND  NOISE 
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SPECTRUM  OF  LOW  FREQUENCY  BROADBAND  NOISE  AS  A FUNCTION  OF  NON  DIMENSIONAL  FREQUENCY 


FIGURE  5.2  THE  EFFECT  OF  LEADING  EDGE  SERRATION  ON  6%  BANDWIDTH  HIGH 
FREQUENCY  BROADBAND  NOISE  SPECTRA 


LENGTH  OF  SERRATION  (%  OF  SPAN) 
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SOUND  PRESSURE  LEVEL,  db  ref.  0.0002 


ROTOR  THRUST  #- 


FIGURE  5.6  THE  EFFECT  OF  ROTOR  COLLECTIVE  PITCH  ON  6 % BANDWIDTH 
HIGH  FREQUENCY  BROADBAND  NOISE  SPECTRA. 


SOUND  PRESSURE  LEVEL  db,  Ref.  0.0002  /xbars 


SOUND  PRESSURE  LEVEL  d b,  Ref.  0.0002  h bars 
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FREQUENCY,  kHZ 

FIGURE  5.8  EFFECT  OF  ADVANCE  RATIO  ON  HIGH  FREQUENCY 
BROADBAND  NOISE. 


FIGURE  5.9  VARIATION  OF  PEAK  INTENSITY  OF  HIGH  FREQUENCY 
BROADBAND  NOISE  WITH  ADVANCE  RATIO. 
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FIGURE  5.11  VARIATION  OF  PEAK  INTENSITY  OF  HIGH  FREQUENCY  BROADBAND  NOISE  WITH 
NUMBER  OF  BLADES. 


FIGURE  5.13  WAKE  FLOW  PATTERN  BEHIND  AN  AIRFOIL. 


Frequency  in  kHz 


FIGURE  5.16 

TYPICAL  SPECTRA  OF  HIGH  FREQUENCY  BROADBAND  NOISE  AT 
VARIOUS  ROTATIONAL  SPEEDS 


FREQUENCY,  kHz. 


% TURBULENCE  INTENSITY 
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FIGURE  6.1 

COMPARISON  OF  THE  PREDICTED  PEAK  FREQUENCY  WITH  MEASURED 
FULL  SCALE  ROTOR  DATA 


COMPARISON  BETWEEN  PREDICTED  AND  MEASURED  FULL  SCALE  ROTOR  SOUND  PRESSURE  LEVELS 


